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Introduction 


Microwave  energy  has  the  advantages  that  at  low  power  levels  there  are  no  radiation 
dangers,  no  contrast  agents,  and  the  examinations  are  comfortable.  This  research 
addresses  directly  the  decision-theoretic  task  of  detection  and  localization  of  breast 
tumor,  using  microwave  diffraction  measurements.  Bayesian  detection  theory  is  used  to 
improve  the  probability  of  correct  detection  and  localization.  This  improvement  in 
performance  is  possible  because  conventional  imaging  techniques,  by  themselves,  usually 
emphasize  resolution  and  contrast,  and  leave  the  incorporation  of  uncertainties  and 
decisions  primarily  to  algorithms  or  human  observers  that  post  process  the  reconstructed 
image.  This  approach  augments  conventional  medical  image  processing  and  provides 
additional  processing  of  the  scattered  microwave  field  to  aid  the  radiologist  in  dealing 
with  uncertainties  that  are  an  inherent  part  of  the  decisions. 

Sha  et  al  have  compiled  results  from  experiments  in  the  literature  that  show  that  the 
microwave  dielectric  properties  of  malignant  tissue  are  different  from  those  of  normal 
breast  tissue.  Signal  detection  theory,  in  its  most  fundamental  form,  provides  a 
framework  for  incorporating  this  knowledge  of  breast  tissue  characteristics  directly  into 
the  design  of  optimal  task  oriented  information  processing  algorithms  to  aid  in  the 
detection  of  breast  cancer.  In  addition,  signal  detection  theory  enables  one  to  obtain  upper 
limits  of  detection  and  localization  performance  as  a  function  of  uncertainties  in  the 
microwave  properties  of  breast  tissues,  using  quantitative  measures  such  as  the  ROC 
(probability  of  detection  vs.  false  alarm)  and  PCL  (probability  of  correct  localization). 

The  proposed  algorithm  uses  the  direct  microwave  diffraction  measurements  and 
incorporates  the  fact  that  spatially  adjacent  tissues  are  similar  in  their  permittivity  values, 
but  normal  and  malignant  breast  tissues  have  high  contrast  in  the  mean.  These  research 
results  demonstrate  the  advantages  of  incorporating  the  microwave  diffraction 
measurements  directly  into  the  computer-aided  algorithm  design. 

Body 

Detailed  descriptions  of  the  research  done  on  this  project  during  this  reporting  period  are 
summarized  below  in  References  [1-4].  Copies  of  the  conference  papers  and  posters 
cited  [1-4]  are  included  in  the  appendix  of  this  report. 

Background  —  Most  of  the  past  research  in  this  field  has  focused  on  the  study  of  the 
dielectric  properties,  the  design  of  the  microwave  imaging  prototypes,  and  the 
improvement  of  the  EM  forward  and  inverse  algorithms.  However,  none  of  this  research 
has  incorporated  signal  detection  theory  directly  into  the  microwave  imaging  at  the 
measurement  level.  Markov  Random  Fields  (MRF)  and  detection  theory  have  been 
applied  in  mammography  for  diagnosis.  However  these  only  assumed  a  simple 
deterministic  disk  object  model  or  did  not  utilize  the  a  priori  knowledge  of  the 
projections.  This  research  presents  Bayesian  algorithms  for  Scattered  Electromagnetic 
fields  through  an  Uncertain  Permittivity  Image  which  incorporates  the  knowledge  of  the 
a  prior  permittivity  image  modeled  by  the  MRF,  the  measurement  noise,  as  well  as  the 
physical  model  of  the  forward  scattered  electric  field.  The  Bayesian  algorithms  for  the 
Uncertain  Permittivity  Image  and  the  Threshold  Image  Processors  are  also  presented  for 
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comparisons.  Our  approach  exploits  the  propagation  of  the  scattered  microwave  fields  to 
develop  better  diagnostic  decision  aids.  Additional  references  to  this  background 
information  are  contained  in  Ref  3 

On  this  concept  research  award,  Sha,  the  graduate  student  on  this  project,  with  the  help  of 
two  other  students,  has  compiled  results  from  experiments  in  the  literature  that  show  that 
the  microwave  dielectric  properties  of  malignant  tissue  is  different  from  that  of  normal 
breast  tissue  [2],  Hence  microwave  imaging  has  the  potential  of  providing  a  tool  for 
improving  the  diagnosis  of  breast  cancer  with  no  known  radiation  dangers  [2]. 

Although  the  experimental  literature  has  provided  promising  evidence  that  the  dielectric 
properties  of  malignant  and  benign  tissue  are  different,  this  project  investigated  the 
impact  of  this  on  the  ability  to  detect  and  localize  tumors.  Since  tissue  properties  change 
spatially,  statistical  models  were  developed  on  this  project  that  model  simple 
uncertainties  in  the  tissue  permittivity  within  the  framework  of  microwave  imaging[3,  4], 
The  Markov  Random  Field  was  used  to  model  the  breast  permittivity  cross  section  as  a 
propagating  medium,  and  this  was  incorporated  into  the  forward  Electromagnetic  (EM) 
propagation  model  [1],  to  predict  the  random  field  of  the  EM  measurements  at  a  received 
array  of  sensors.  These  statistical  models  of  tissue  permittivity  were  then  incorporated 
into  an  optimal  signal  detection  theory  framework  in  which  task-oriented  goals  such  as 
detection  and  localization  of  a  tumor  drive  the  information  processing  [3,  4],  Using 
Bayesian  signal  detection  theory,  the  likelihood  ratio  for  tumor  detection  and  the  a 
posteriori  probability  of  tumor  location  were  computed  and  displayed.  The  ROC 
(receiver  operating  characteristics)  was  used  as  a  quantitative  performance  measure. 

These  results  provide  an  upper  bound  on  the  task-driven  goal,  namely  the  detection  of  the 
presence  or  absence  of  a  tumor  as  a  function  of  tumor  permittivity  contrast,  size,  noise, 
and  local  spatial  permittivity  uncertainties  that  characterize  the  tissue  in  microwave 
imaging  [3,  4],  In  addition,  the  probability  of  correct  localization  (PCL)  was  obtained  as 
a  quantitative  measure  of  how  well  one  can  determine  the  location  of  a  tumor,  as  a 
function  of  tumor  contrast,  size,  and  spatial  tissue  properties  [3,4], 

Details  of  this  research  supported  by  this  project  are  presented  in  the  conference  papers 
and  poster  sessions,  references  [1-3],  with  copies  of  drafts  of  these  references  [1-3] 
included  in  the  appendix. 


Key  Research  Accomplishments 

•  Experimental  results  were  compiled  from  the  literature.  Those  results  show  that  the 
dielectric  properties  of  malignant  and  normal  breast  tissue  at  microwave  frequencies 
are  different.  This  compilation  was  led  by  Liewei  Sha,  a  graduate  student  on  this 
project,  and  the  results  are  published  in  Reference  [2],  Liewei  Sha,  Erika  Ward ,  and 
Brandon  Story  ,"A  review  of  dielectric  properties  of  normal  and  malignant  breast 
tissue,"  Proceedings  of  the  IEEE  SoutheastCon  2002,  pp.  457  -462,  Columbia,  South 
Carolina,  April  5-7,  2002. 
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•  Developed  Markov  Random  field  statistical  models  that  characterize  the 
uncertainties  in  the  spatial  properties  of  tissue  permittivity  at  microwave  frequencies. 
Incorporated  this  information  into  the  forward  Electromagnetic  (EM)  propagation  to 
predict  the  random  field  of  the  EM  measurements  at  a  received  array  of  sensors  [3,4]. 

•  Used  an  optimal  signal  detection  theory  framework  in  which  task-oriented  goals  such 
as  detection  and  localization  of  a  tumor  drive  the  information  processing.  Bayesian 
approaches  were  developed  to  compute  the  likelihood  ratio  for  tumor  detection  and 
the  a  posteriori  probability  display  of  tumor  location.  [3,4] 

•  Although  the  permittivity  between  malignant  and  benign  tissue  has  been  reported  in 
the  experimental  literature  to  be  different,  the  question  is  how  does  this  translate  to 
ones  ability  to  detect  malignant  tissue.  Using  the  optimum  signal  detection  and 
localization  information  processing  algorithms  derived  above,  ROC’s  were  obtained 
to  provide  an  upper  bound  on  the  detection  of  the  presence  or  absence  of  a  tumor. 
Results  were  obtained  as  a  function  of  tumor  permittivity  contrast  between  malignant 
and  normal  tissue,  size,  noise,  and  local  spatial  permittivity  uncertainties  that 
characterize  the  tissue  at  microwave  frequencies.  [3,4] 

•  The  probability  of  correct  localization,  PCL,  was  also  obtained  as  a  quantitative 
measure  of  how  well  one  can  determine  the  location  of  a  tumor,  as  a  function  of 
tumor  contrast,  size,  tissue  spatial  permittivity  characteristics,  and  noise.  [3,4] 

Reportable  Outcomes 


Liewei  Sha,  Erika  Ward ,  and  Brandon  Story  ,"A  review  of  dielectric  properties  of 
normal  and  malignant  breast  tissue,"  Proceedings  of  the  IEEE  SoutheastCon  2002,  pp. 
457  -462,  Columbia,  South  Carolina,  April  5-7,  2002. 

Liewei  Sha,  Loren  W.  Nolte,  Zhong  Qing  Zhang  and  Qing  H.  Liu,  "Performance  analysis 
for  Bayesian  microwave  imaging  in  decision  aided  breast  tumor  diagnosis,”  Proceedings 
of  the  2002  TEF.F.  International  Symposium  on  Biomedical  Imaging  Washington  DC., 
pp.  1039-1042,  Washington,  DC,  July7-10,  2002. 


Liewei  Sha,  Loren  W.  Nolte,  Zhong  Qing  Zhang  and  Qing  Huo  Liu,  "Decision 
Aided  Algorithms  for  Breast  Tumor  Diagnosis  Using  Microwave  Diffraction 
Measurements,"  Poster,  Era  of  Hope  2002  Department  of  Defense  Breast  Cancer 
Research  Meeting,  Orlando,  Florida,  September  25-28,  2002. 


Conclusions 


Decision-aided  breast  tumor  diagnostic  algorithms  based  on  signal  detection  theory  and 
microwave  energy,  which  has  no  radiation  danger  at  low  power  levels,  have  the  potential 
of  providing  additional  information  for  radiologist  so  as  to  improve  the  probability  of 
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early  detection  of  breast  tumors  as  well  as  their  correct  locations.  This  Bayesian 
microwave  imaging  approach  is  driven  by  the  ultimate  decision  tasks  of  whether  or  not  a 
tumor  is  present,  and  if  so  where  is  it.  It  benefits  from  incorporating  the  a  priori 
knowledge,  although  uncertain,  of  the  permittivity  properties  of  normal  and  malignant 
breast  tissue.  This  approach  also  gains  a  decision  performance  advantage  by  processing 
the  measurements  directly  using  signal  detection  theory,  rather  than  post  processing  a 
reconstructed  image  where  some  of  the  information  needed  for  the  decision  has  been  lost. 
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Active  Microwave  Imaging  I — 2-D  Forward  and 
Inverse  Scattering  Methods 

Qins  Huo  Liu,  Senior  Member,  IEEE,  Zhong  Qing  Zhang,  Senior  Member,  IEEE,  Tonghui  T.  Wang, 

J  A.  Bryan,  Gary  A.  Ybarra,  Member,  IEEE,  Loren  W.  Nolte,  Life  Senior  Member,  IEEE,  and 
William  T.  Joines,  Life  Senior  Member,  IEEE 


Abstract — Active  microwave  imaging  (MWI)  for  the  detection 
of  breast  tumors  is  an  emerging  technique  to  complement  existing 
X-ray  mammography.  The  potential  advantages  of  MWI  arise 
mainly  from  the  high  contrast  of  electrical  properties  between 
tumors  and  normal  breast  tissue.  However,  this  high  contrast  also 
increases  the  difficulty  of  forming  an  accurate  image  because  of 
increased  multiple  scattering.  To  address  this  issue,  we  develop  fast 
forward  methods  based  on  the  combination  of  the  extended  Born 
approximation,  conjugate-  and  biconjugate-gradient  methods, 
and  the  fast  Fourier  transform.  We  propose  two  nonlinear  MWI 
algorithms  to  improve  the  resolution  for  the  high-contrast  media 
encountered  in  microwave  breast-tumor  detection.  Numerical 
results  show  that  our  algorithms  can  accurately  model  and  invert 
for  the  high-contrast  media  in  breast  tissue.  The  outcome  of  the 
Inversion  algorithms  is  a  high-resolution  digital  image  containing 
the  physical  properties  of  the  tissue  and  potential  tumors. 

Index  Terms — Biomedical  applications,  fast  algorithms, 
imaging,  inverse  scattering,  microwave  techniques. 

I.  Introduction 

EARLY  treatment  of  breast  cancer  can  be  highly  effec¬ 
tive  and  can  significantly  increases  long-term  survival 
[1].  Among  current  clinical  and  experimental  breast-cancer 
imaging  technologies  (see,  e.g.,  [2]),  the  most  popular  method 
for  breast  cancer  screening  and  diagnosis  is  based  on  X-ray 
mammography.  In  spite  of  its  high  resolution,  however,  X-ray 
mammography  has  the  following  shortcomings  [3],  [4]. 

1)  It  has  difficulties  detecting  breast  tumors  at  their  earlier 
stages. 

2)  It  has  a  decreased  effectiveness  in  cases  of  women  with 
dense  breasts  and  has  difficulties  detecting  tumors  located 
near  the  chest  wall  or  underarm. 

3)  It  has  a  limited  specificity  between  benign  and  malignant 
tissue  resulting  in  high  false-alarm  rates  at  clinically  ac¬ 
cepted  detection  rates. 

4)  It  uses  ionizing  radiation. 

5)  There  is  a  discomfort  to  patients  because  of  breast  com¬ 
pression. 
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As  an  alternative  and  complementary  modality  for  breast 
imaging,  microwave  (MW)  techniques  have  been  proposed 
over  the  last  few  years  [5]-[8]  because  of  its  potentially  high 
specificity  for  breast  cancer  diagnosis  due  to  the  high  contrast 
in  electrical  properties  between  normal  and  malignant  human 
breast  tissues. 

Due  to  significantly  different  sodium  concentrations,  fluid 
contents,  and  electrochemical  properties,  a  significant  contrast 
exists  in  electrical  properties  at  MW  frequencies  between 
the  normal  and  malignant  human  breast  tissues  [9]— [13].  For 
example,  at  800  MHz,  the  relative  permittivity  (er)  and  elec¬ 
trical  conductivity  (cr)  for  normal  mammary  tissues  are  around 
er  «16  and  a  «0.16  S/m,  respectively,  while  they  are  er  «57.2 
and  <7  «1.08  S/m,  respectively,  for  a  malignant  breast  tumor. 
The  contrast  is  3.75  for  the  relative  permittivity  and  6.75  for  the 
electrical  conductivity.  This  high  contrast  gives  rise  to  a  large 
electromagnetic  scattering  signal  when  electromagnetic  waves 
are  applied  to  a  malignant  tumor  embedded  in  a  normal  tissue. 

Compared  to  X-ray  mammography,  MW  mammography  has 
the  following  potential  advantages. 

1)  The  high  contrast  in  electrical  properties  can  potentially 
provide  a  high  specificity  to  distinguish  between  normal 
and  malignant  breast  tissues.  MW  mammography  may  be 
able  to  detect,  at  an  earlier  stage,  small  tumors  otherwise 
undetectable  by  X-ray  mammography. 

2)  Microwave  imaging  (MWI)  utilizes  nonionizing  radia¬ 
tion  for  breast  cancer  imaging.  The  MW  radiation  levels 
required  are  well  within  accepted  safety  standards,  and 
there  are  no  contrast  agents  involved. 

3)  The  examinations  using  MWs  are  well  tolerated  since  no 
breast  compression  is  required  for  the  patient  The  MWI 
method  provides  a  useful  alternative  and  complement  to 
the  diagnostic  information  provided  by  X-ray  mammog¬ 
raphy. 

During  the  last  decade,  there  has  been  significant  progress 
in  the  development  of  experimental  prototypes  of  two-dimen¬ 
sional  (2-D)  and  three-dimensional  (3-D)  MWI  systems  based 
on  various  ideas  [5]-|7],  ll5l  systems  image  2-D 
(and  more  recently,  3-D)  small-  and  medium-size  objects,  and 
sometimes  even  the  full  body  [16]— [18].  Substantial  research 
has  also  been  carried  out  specifically  for  breast  cancer  diag¬ 
nosis  with  narrow-band  [14]  and  wide-band  [5],  [6]  signals.  A 
recent  report,  apparently  the  first  clinical  experiments,  shows 
very  promising  results  for  MW  breast  imaging  [7]. 

From  the  practical  viewpoint,  one  of  the  most  significant 
problems  for  current  MWI  systems  is  their  lower  spatial 
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Fig.  1.  Typical  geometry  for  forward  and  inverse  problems  in  MW  breast 
imaging.  Transmitting  and  receiving  antennas  are  mounted  on  the  rim  of  a 
fluid-holding  box  with  an  absorptive  backing. 


resolution*  There  has  been  great  interest  in  inverse  problems 
for  biomedical  applications  of  MWs  to  improve  the  resolution 
[19],  [22],  [23],  [20],  [21].  Unlike  X-ray  imaging  techniques, 
which  only  need  to  consider  the  attenuation  of  the  X-ray,  for 
MWI  technology,  we  need  to  consider  multiple  scattering  (i.e., 
diffraction)  of  MWs  within  the  tissues.  This  is  a  significant 
effect  and  cannot  be  neglected  [18].  From  this  viewpoint, 
the  very  advantage  of  high  contrast  in  electrical  properties 
between  normal  and  malignant  tissues  can  also  pose  a  major 
limitation  if  not  treated  correctly.  This  is  because  the  high 
contrast  introduces  large  multiple  scattering  effects,  leading  to 
a  low  resolution  if  linear  or  quasi-linear  methods  are  used  to 
invert  for  the  object  properties. 

In  this  paper,  we  develop  fast-forward  and  inverse-scattering 
algorithms  to  fully  unravel  multiple  scattering  effects.  These  al¬ 
gorithms  are  applied  to  form  high-resolution  images  for  an  ex¬ 
perimental  prototype  of  a  MW  breast-imaging  system  recently 
developed  at  Duke  University,  Durham,  NC.  Our  imaging  tech¬ 
nique  uses  the  fast  Fourier  transform  (FFT),  the  extended  Bom 
approximation  (EBA),  conjugate  and  biconjugate  gradients,  and 
contrast  source-inversion  (CSI)  methods. 

This  paper  is  organized  as  follows.  We  briefly  present  the  in¬ 
tegral  equation  for  the  2-D  TM  problem  in  Section  II.  The  for¬ 
ward  and  inverse  solution  methods  are  presented  in  Sections  III 
and  IV.  Numerical  results  are  shown  in  Section  V  to  demonstrate 
the  efficacy  of  the  technique.  Conclusions  are  presented  in  Sec¬ 
tion  VI. 


n.  Integral  Equation 

In  this  paper,  we  focus  on  2-D  MWI  problems.  We  assume  a 
TM*  wave  from  a  finite  source  impinging  on  an  inhomogeneous 
medium  with  a  high  contrast  in  both  dielectric  constant  and  con¬ 
ductivity  with  respect  to  a  homogeneous  background  medium. 
In  MW  breast  imaging,  a  matching  material  with  a  nonreflec- 
tive  backing  is  usually  used  to  eliminate  the  reflections  at  the 
tissue/air  interface,  as  illustrated  in  Fig.  1 .  In  that  case,  the  back¬ 
ground  medium  is  the  matching  fluid  with  properties  similar  to 
a  normal  breast  tissue. 

Fig.  1  illustrates  the  general  scenario  of  the  problem,  where 
an  inhomogeneous  object  with  dielectric  constant  tr(r),  con¬ 
ductivity  a( r),  and  a  constant  magnetic  permeability  fib  is  im¬ 


mersed  in  a  background  medium  with  the  corresponding  pa¬ 
rameters  cr&,  <r$,  and  fib-  With  an  implied  time  dependence  of 
exp  (jut),  the  complex  permittivity  for  the  object  and  back¬ 
ground  can  be  written  as  e  =  eo€r  —  jcrju  and  tb  =  eo^rfc  ~ 
job/ o>.  In  the  2-D  TM*  case,  all  field  and  source  variables  are 
only  functions  of  r  =  ( x ,  y ).  For  an  electric  current  source 
J  ss  zJz(r)  in  such  an  inhomogeneous  medium,  the  electric 
field  E  =  zEz{ r)  satisfies  the  integral  equation 

Ez{r)  =  ETi r)+ Jdr'g(T-T')  [**(r')-*?]l?.(*')  (D 

D 

where  the  complex  wavenumbers  k  and  kb  are  defined 
by  k2  =  u2fibt  and  fcj  =  respectively,  and 

g  (r  -  r')  =  (l/4j)ifo2)  {h\r  -  r'|)  is  the  Green’s  function 
for  a  homogeneous  background  medium.  The  incident  field 
Efc  in  (1)  is  the  field  in  the  background  medium,  i.e.,  in  the 
absence  of  the  scatterer  (or  k  =  kb). 

Equation  (1)  is  the  central  equation  we  use  in  this  paper  for 
both  the  forward  and  inverse  solutions.  In  the  forward  problem, 
both  ffie  medium  properties  and  the  domain  of  inhomogeneity 
D  are  known;  our  aim  is  to  solve  (1)  to  obtain  the  electric  field. 
In  the  inverse  problem,  the  electromagnetic  fields  are  measured 
at  some  discrete  points,  usually  outside  the  domain  of  inhomo¬ 
geneity;  the  medium  properties  (including  the  domain  of  inho¬ 
mogeneity)  are  unknowns  to  be  determined.  However,  we  as¬ 
sume  that  the  scatterer  is  finite  and,  thus,  can  always  be  enclosed 
in  a  large  rectangular  region  D. 

A.  Historical  Note 

There  is  a  large  body  of  literatures  dealing  with  both  the  for¬ 
ward  and  inverse  solutions  for  integral  equations  similar  to  (1). 
It  is  not  our  purpose  to  present  a  thorough  review  of  these  works. 
Here,  we  will  briefly  summarize  the  developments  leading  to  the 
methods  used  in  this  study. 

In  the  context  of  low-frequency  electromagnetic  induction 
applications,  the  forward  and  inverse  problems  were  solved 
with  the  EBA  in  [25]-[27].  An  FFT  accelerated  EBA  method 
was  developed  as  fast-forward  and  inverse  solvers  in  [28],  and 
has  been  used  as  a  preconditioner  for  the  conjugate-gradient 
fast  Fourier  transform  (CG-FFT)  method  in  [29],  [30],  and 
[45].  However,  the  application  in  MWI  involves  significantly 
different  physics  than  in  our  previous  low-frequency  problems 
where  the  wave  fields  are  diffusive  because  of  the  dominant 
conduction  current  The  dominant  wave  phenomenon  in  the 
MWI  application  makes  both  the  forward  and  inverse  problems 
more  challenging.  The  biconjugate-gradient  fast  Fourier  trans¬ 
form  (BCG-FFT)  method  used  in  this  paper  is  an  application 
of  our  recent  method  presented  in  [49]. 

HI.  Forward-Solution  Methods 

There  are  two  purposes  for  the  forward  solution.  First,  it  will 
be  used  to  simulate  and  calibrate  the  imaging  system.  Secondly, 
it  will  be  used  to  test  the  inverse  algorithms  by  providing  the 
simulated  “measured”  data.  Furthermore,  other  inverse  algo¬ 
rithms  (unlike  those  in  this  paper)  may  require  repeated  forward 
solutions. 
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In  the  forward  problem,  the  unknown  electric  field  Ez  ap¬ 
pears  both  on  the  left-hand  side  and  inside  the  integrand  of  (1). 
The  conventional  way  of  solving  (1)  is  the  method  of  moments 
(MOM),  which  discretizes  the  integral  equation  with  N  electric 
field  unknowns  in  terms  of  basis  functions  inside  D ,  and  solves 
the  resulting  matrix  equation  by  direct  matrix  inversion.  Unfor¬ 
tunately,  this  direct  method  is  prohibitively  expensive,  costing 
0(iV3)  CPU  time  and  0(N2)  computer  memory. 

A .  EBA 

An  efficient  way  to  solve  (1)  is  to  use  the  Bom  approximation 

[24] ,  i.e.,  to  assume  that  the  electric  field  inside  the  integrand  in 
(1)  can  be  approximated  by  the  incident  field  E^c.  The  total 
field  can  then  be  easily  calculated  by  performing  the  integral  of 
the  Green’s  function  operating  on  the  Bom  approximated  cur¬ 
rent  density. 

The  Bom  approximation  has  a  very  limited  range  of  validity 
when  the  contrast  and/or  the  size  of  the  scatterer  is  small 

[25] — [27]-  To  increase  the  range  of  validity,  one  can  use  the 
EBA  put  forward  in  [25].  In  essence,  the  EBA  makes  use  of 
the  fact  that  the  Green’s  function  g  (r  -  r')  is  highly  peaked 
(in  feet,  singular)  when  r  — >  r',  and  becomes  relatively  small 
for  larger  |r  -  r'|.  Hence,  one  can  replace  Ez  (r')  inside  the 
integrand  in  (1)  with  Ez{t)  and  approximate  (1)  by 

Ejft r)  »  £“c(r)+££(r)A|  J dr'g( r-r')x(r'),  t  €  D 

°  (2) 
where  x  =  (k2/kl  ~  l)  is  the  contrast  function.  Consequently, 
the  electric  field  in  the  integrand  of  (1)  can  be  more  accurately 
approximated  by  the  solution  of  (2)  when  this  approximated 
field  JEf  from  (2)  is  substituted  into  the  integrand  of  (1)  to  yield 

£?D)(r) 

=  £T(r)  +  kl  J dr'g(r  -  r')x(r') M"1  (r')2?r(r') 

D 

(3) 

where 

M(r)  =  1  -  k\  J  dr'g(r  -  r')x(r') ,  r  €  D.  (4) 

D 

It  has  been  shown  that  the  EBA  has  a  significantly  larger 
range  of  applicability  than  fee  Bom  approximation  [25],  [26], 
[28],  [29],  [30].  Furthermore,  the  computational  cost  for  (3) 
remains  essentially  the  same  as  the  Bom  approximation,  i.e., 
O(MfN)  arithmetic  operations,  where  Mf  is  the  number  of 
field  points.  However,  (4)  requires  0(N2)  arithmetic  operations 
since  M(r)  is  needed  for  all  N  points  within  the  object.  Thus, 
the  total  cost  of  the  EBA  is  0{C\MfN  +  C^iV2),  where  C\ 
and  C2  are  constants.  This  is  still  much  more  efficient  than  the 
direct  MOM  with  0(Ar3)  arithmetic  operations. 

B.  Using  FFT  to  Speed  up  the  EBA :  FFT-EBA  Method 

We  adopt  an  efficient  method  for  the  electromagnetic  induc¬ 
tion  problem  [28]— [30]  to  speed  up  the  EBA  computation  of  the 
electric  field  at  Ar  points  within  the  scatterer.  This  method  is 


based  on  the  observation  that  the  integrals  in  (3)  and  (4)  are  con¬ 
volutional.  Hence,  by  using  the  convolution  theorem,  one  can 
solve  for  the  electric  field  within  the  EBA  approach  as 

£fD)(  r) 

(5) 

where  F  and  T~l  stand  for  the  2-D  forward  and  inverse 
Fourier  transforms,  respectively,  which,  in  their  discrete  forms 
with  zero-padded  arrays,  can  be  achieved  efficiently  through 
the  FFT  algorithms  with  0(Arlog2  Ar)  arithmetic  operations. 
Equation  (5)  also  defines  the  operator  £eb- 

This  improvement  reduces  the  computational  cost  of  the 
EBA  from  0(N2)  to  0(N  log2  AT).  Furthermore,  it  makes  pos¬ 
sible  the  following  hybridization  of  the  EBA  with  the  CG-FFT 
and  BCG-FFT  methods  without  increasing  the  overall  cost  to 
0(N2).  For  the  convenience  of  the  following  discussions,  we 
refer  to  this  improved  approach  as  the  FFT-EBA  method. 

C.  EBA  as  a  Preconditioner  for  CG—FFT  and  BCG— FFT 
Methods 

Although  the  EBA  is  more  accurate  than  the  Bom  approxi¬ 
mation,  its  error  increases  when  the  object  becomes  very  large 
in  size  or  in  contrast  against  the  background.  An  efficient 
full-wave  method  that,  in  principle,  does  not  suffer  from  this 
limitation  is  the  CG-FFT  method  [33]— [36].  This  iterative 
method  is  more  efficient  than  the  MOM,  as  it  requires  only 
0(K N  log2  N)  arithmetic  operations,  where  K  is  the  number 
of  CG  iterations.  In  our  previous  studies  of  the  electromagnetic 
induction  problem  [29],  [30],  we  propose  to  incorporate  the 
EBA  to  further  improve  the  efficiency  of  the  CG-FFT  or 
BCG-FFT  methods  [49]. 

Our  aim  here  is  to  solve  (1)  for  the  unknown  electric  field  Ez 
within  the  inhomogeneous  region.  Noting  that  the  integral  in  (1) 
is  a  2-D  convolution  between  the  Green’s  function  #(r)  and  in¬ 
duced  current  source  x(r)2?z(r),  we  can  invoke  the  convolution 
theorem  to  rewrite  (1)  in  an  operator  form 

C{EZ]  =  Ez  -  klF~l  {jf  [<?(r)]  T  [x{r)Ez] }  =  £fc  •  (6) 

This  equation  can  then  be  solved  iteratively  by  the  CG-FFT 
method  [33],  [34],  [36]. 

We  apply  the  EBA  as  the  preconditioner.  Hence,  instead  of 
solving  (6)  directly,  we  solve  an  equivalent  problem 

C^CEZ  =  (7) 

where  the  preconditioner  is  the  EBA  operator  defined  in  (5). 
Since  the  preconditioner  is  achieved  by  the  FFT  algorithm 
in  (5),  this  new  (7)  can  again  be  solved  efficiently  by  the 
CG-FFT  method.  With  this  preconditioning,  this  scheme  con¬ 
verges  much  fester  than  the  regular  CG—FFT  method  [28],  [29]. 
In  reality,  we  found  that  even  using  the  EBA  preconditioner 
in  only  the  first  step  of  the  conjugate-gradient  (CG)  method 
can  significantly  improve  the  convergence  speed  for  typical 
electromagnetic  problems.  This  partial  preconditioning  avoids 
the  additional  cost  of  the  preconditioning  in  the  majority  of  the 
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CG  steps.  The  following  numerical  results  are  obtained  by  this 
preconditioning  scheme. 

Equation  (7)  can  also  be  solved  by  using  the  biconjugate-gra¬ 
dient  (BCG)  method  [31],  [32],  [49]  with  the  FFT  algorithm. 
The  resulting  EBA  preconditioned  BCG-FFT  method  is  more 
efficient  than  the  CG-FFT  method  described  above.  For  proce¬ 
dures  of  the  BCG  iterations,  see,  for  example,  [31]  and  [49]. 

Note  that  there  is  no  approximation  in  the  above  EBA 
preconditioned  CG-FFT  and  BCG-FFT  methods  aside  from 
the  discretization  procedure.  The  combination  of  the  EBA  and 
CG-FFT  and  BCG-FFT  methods  provides  a  seamless  method, 
which  is  both  efficient  and  accurate  for  low-  and  high-contrast 


scatterers.  In  practice,  we  define  an  acceptable  error  criterion 
emin  based  on  the  i/2-norm  residual  error 
|| £EZ  -  Eh 


C2 


II  & 


— 


to  terminate  the  iteration  process.  For  lower  contrast  problems, 
if  the  initial  step  of  the  EBA  solution  suffices,  no  further  CG  or 
BCG  procedures  are  needed;  for  higher  contrast  problems,  the 
CG  or  BCG  iterations  will  continue  until  the  residual  error  is 
smaller  that  emin.  Note  that  e2  can  always  be  calculated  even 
though  the  exact  solution  is  not  known.  Compared  with  the  reg¬ 
ular  CG-FFT  and  BCG-FFT  methods,  these  EBA  precondi¬ 
tioned  methods  are  more  efficient  [29]. 


this  FFT  enhanced  version  adapted  to  the  current  problem  in 
Cartesian  coordinates  [30]  as  follows. 

Step  1)  The  induced  source:  The  first  step  of  this  scheme  is 
to  invert  for  the  induced  source  function  w(t)  from 
(8).  With  the  trapezoidal  rule,  (8)  can  be  discretized 
into  a  linear  system  of  equations 

A  *  w  =  d  (10) 

where  A  is  an  M  x  N  matrix,  d  is  the  data 
vector  containing  the  measured  scattered  field, 
N  is  the  number  of  unknowns  in  domain  D,  and 
M  —  MrMt  is  the  number  of  measured  data  points 
determined  by  the  number  of  source  locations  Mt 
and  the  number  of  receiver  locations  Mr.  Note  that 
this  equation  for  the  unknown  w  is  a  classical  linear 
ill-posed  problem.  A  minimum-norm  solution  to 
(10)  can  be  obtained  by  an  optimization  procedure. 
That  is,  we  find  a  solution  w  that  will  minimize  the 
1/2  norm  functional 

1  =  ||d-  A*w||2+7lM|2  (11) 

where  7  is  a  regularization  parameter.  Minimizing  I 
yields  a  linear  equation  as  follows: 


IV.  Nonlinear  Inverse-Scattering  Algorithms 

In  the  inverse  problem,  the  unknown  material  properties  er(r) 
and  <r(r)  are  inferred  from  the  scattered  electric  field  measured 
at  some  discrete  locations.  This  inverse  problem  is  nonlinear,  as 
the  material  properties  appear  both  in  the  contrast  function  x 
and  in  Ez  in  (1).  Furthermore,  this  inverse  problem  is  ill  posed 
because  of  the  sparsity  of  measured  data  contaminated  by  noise 
[38]. 

Methods  for  solving  the  electromagnetic  inverse-scattering 
problem  include  linear  [37]  and  nonlinear  [26],  [28],  [38], 
[40],  [42],  [43],  [47],  [48]  inverse  methods.  (See  [44]  for  more 
complete  references.)  In  this  study,  we  apply  two  methods  for 
nonlinear  inverse  problems.  The  first  is  an  improved  two-step 
inversion  method  based  on  the  EBA  [25]-[28],  [45].  This  is 
a  compromise  between  the  linear  and  full  nonlinear  inversion 
methods.  The  second  nonlinear  method  is  based  on  the  CSI 
method  [46]— [48]  with  an  improved  initial  solution  through  the 
two-step  nonlinear  inversion. 

A.  Two-Step  Inverse  Method  Based  on  EBA 

1)  FFT  Enhanced  Two-Step  Inversion  Method:  The  EBA  in 
(3)  can  be  rewritten  for  the  scattered  field  as 

E^a(tr,tt)  =  kl  J  dr'  g(rR  -  r')v}(r')E™(r' ,rT)  (8) 

D 

where,  for  clarity,  we  have  explicitly  included  the  transmitter 
and  receiver  locations  tt  and  tr  in  the  electric  field.  In  the 
above,  w(r)  is  defined  as 

w(r)=X(r)M-1(r)  (9) 

where  M  is  given  in  (4). 

The  two-step  linear  inverse  method  has  been  applied  in  [27] 
and  enhanced  by  the  FFT  algorithm  [28].  Here,  we  summarize 


(A*  •  A  -f  7l)w  =  A*  •  d.  (12) 

Equation  (12)  can  be  solved  iteratively  by  the  CG 
procedure  with  the  cost  of  O(MN)  per  iteration. 

Step  2)  The  contrast  function:  Once  w(t)  is  obtained,  our 
second  linear  inversion  step  is  to  invert  for  x(r)  (and, 
thus,  er(r)  and  a(r))  using  (9)  and  (4).  That  is  to 
solve  the  following  equation: 

X(r)  +  tfw(r)  ■  J7-1  {^T?JF[x(r)] }  =  w(r)-  (l3) 

Equation  (13)  is  a  linear  equation  for  the  unknown 
contrast  function  x(r),  a  function  representing  the 
difference  between  the  anomaly  and  background. 
Equation  (13)  is  a  well-posed  linear  problem.  It  can 
be  solved  iteratively  by  using  the  CG  method  to  find 
unknown  complex  permittivity  c  =  eocr  —  jcr/cj.  In 
the  CG  procedure,  the  Green’s  function  operations 
are  accelerated  by  the  FFT  algorithm.  This  repre¬ 
sents  an  acceleration  over  the  earlier  versions  of  the 
two-step  inversion  scheme  [26],  [27],  and  reduces 
the  number  of  arithmetic  operations  to  0(N  log2  N) 
in  each  CG  iteration  [28],  [30]. 

With  these  procedures,  the  original  nonlinear  in¬ 
verse  problem  has  been  converted  into  two  linear  in¬ 
verse  problems  based  on  the  EBA.  This  scheme  uses 
the  FFT  algorithm  to  accelerate  the  second  linear 
inversion  step  to  greatly  reduce  the  computational 
cost  from  0(jY3)  to  0(N  log2  N),  where  N  is  the 
number  of  unknown  pixels  in  the  inverse  model. 
Next,  we  will  improve  the  memory  efficiency  in  the 
first  linear  inversion  step. 
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2)  Improved  Two-Step  Inversion  Scheme:  The  above 
two-step  inversion  scheme  is  a  significantly  enhanced  ver¬ 
sion  of  earlier  work  since  the  total  CPU  time  is  reduced  to 
0{CxKiMN  +  C2K2N\q%2N),  where  Cx  and  C2  are 
constants,  and  K\  and  K2  are  the  numbers  of  the  CG  iterations 
in  the  two  steps,  respectively  [28].  However,  the  memory  re¬ 
quirement  is  still  large  for  the  first  step  as  A  requires  0(M N ) 
storage.  We  can  further  improve  the  memory  requirement  in 
a  way  similar  to  the  CSI  method  [47],  [48].  Defining  S  and 
D  as  the  domain  of  measured  data  and  the  imaged  domain, 
respectively,  and  Green’s  operator  Gs  and  Gjj  as 

Gs,d[-]  =kb  J  dr'9  (r  “  r')  H)  r  e  5  or  D  (14) 

D 

we  can  rewrite  (8)  as 

E?  (rfi!rr)  =  {Gs[£“c(rVTMr')]  }(r«).  (15) 

The  first  step  of  the  EBA  inversion  for  ti/(r)  can  then  be  stated 
as  the  minimization  of  the  functional 

F  =  || Ef  -  GS [E?cw]  II’  +  7|MfD  (16> 

where  S  is  the  measured  data  domain.  This  functional  is  basi¬ 
cally  the  same  as  that  in  (1 1),  except  that  it  is  in  an  undiscretized 
form. 

Since  (16)  is  a  quadratic  functional  of  w,  it  can  be  minimized 
by  the  conjugate  gradient  procedure.  Defining  the  data  error 

Pj  =  Elf  -  GS  [E™w] 

we  can  update  wn  by  the  following  CG  procedure: 

Wn  =  ltfn-1  +  CtnVn  (17) 


where  the  update  directions  { vn }  are  chosen  by  the 
Polak-Ribiere  CG  directions 


vq  =0 


vn  —dn  + 


(dn,dn  —  l)l> 

{dn-udn-l  )D  " 


(18) 

n  >  1.  (19) 


The  gradient  (Frechet  derivative)  of  the  cost  functional  with  re¬ 
spect  to  w  is  given  by 


the  inversion  is  much  faster  than  the  CSI  because  the  system  of 
equations  for  w  is  much  smaller. 

Note  that  even  though  the  CPU  time  for  this  modified  first 
step  remains  0(AiMiV),the  computer  memory  cost  is  reduced 
from  0(MAr)  to  0(MtN),  where  M  =  MTMR  is  the  number 
of  measured  data.  After  w  is  obtained,  we  solve  (13)  to  deter¬ 
mine  the  complex  permittivity  profile  with  the  computational 
complexity  0(Arlog2  A’). 

B.  EBA  Preconditioned  CSI 

For  many  problems,  the  two-step  inversion  scheme  provides 
satisfactory  results.  However,  since  the  EBA  is  an  approximate 
method,  the  inverse  procedure  based  on  the  EBA  to  reconstruct 
the  complex  permittivity  profile  may  not  always  be  adequate, 
especially  for  very  high  contrasts.  Under  those  circumstances, 
some  nonlinear  inversion  methods  can  be  employed  to  achieve 
better  results.  The  recently  developed  CSI  method  has  an  im¬ 
portant  feature  in  that  it  does  not  require  a  forward  solution  in 
the  inverse  iterations  [47],  [48].  We  adopt  this  CSI  method  for 
MW  breast  imaging  by  defining  an  object  functional 

F  =  vs  £  ||£T*  -  os  M  ||“  +  TO  £  Mi”  -  •> 

3  3 

+xGD[wj\\\2D 

(22) 

where  Wj  =  x(T)E*  (r>  tt.j)  for  r  €  D  is  the  contrast  source, 
and  the  normalization  constants  rjs  and  t)d  are  chosen  so  that 
both  terms  in  (22)  are  equal  to  one  if  wj  ~  0.  The  procedures  to 
minimize  F  in  (22)  for  the  solution  of  x  can  be  found  in  [48]. 
It  is  worth  noting  that  even  though  (22)  is  similar  to  (16),  the 
minimization  of  (22)  is  more  time  consuming  since,  here,  Wj  is 
a  function  of  the  source  location. 

Given  the  measured  scattered  electric  field,  we  can  invert  for 
the  permittivity  profile  by  the  CSI  method.  However,  in  contrast 
to  the  earlier  CSI  implementations,  we  use  the  above  two-step 
inversion  result  rather  than  the  back-propagation  as  the  initial 
solution  to  speed  up  the  convergence.  Preliminary  results  show 
that  this  improved  CSI  procedure  accelerates  the  convergence 
of  the  inversion  as  the  two-step  EBA  inversion  already  achieves 
a  reasonable  approximation.  For  our  typical  examples,  the  EBA 
preconditioned  CSI  method  requires  less  than  half  of  the  itera¬ 
tions  and  CPU  time  of  the  original  CSI  method. 


iWT  r  t 

dn  =  -K2  £  [i?r  (r,  r TJ)]  J  9* (r  -  r')py,„_i (r')dr' 

3=1  JS 

+7™n-l>  T  <=D.  (20) 


where  the  asterisk  denotes  the  complex  conjugate.  By  mini¬ 
mizing  F  in  (16),  the  constant  an  is  determined  as 


(dn,Vn)D 


||Gs[£rt;*]||s+7|MD 


(21) 


Note  that  an  important  difference  between  this  procedure  and 
that  in  the  CSI  (see  [47])  is  that  w  in  (16),  unlike  the  contrast 
source  in  CSI,  is  not  a  function  of  source  location.  Therefore, 


V.  Numerical  Results  and  Discussions 

Our  numerical  models  will  aim  to  simulate  a  prototype  MWI 
system  for  breast  cancer  imaging  developed  at  Duke  Univer¬ 
sity.  Fig.  1  schematically  shows  the  basic  geometry  modeled 
in  this  study.  The  rectangular  cylinder  is  filled  with  a  matching 
fluid  whose  electrical  properties  are  close  to  those  of  the  normal 
breast  tissue.  The  outer  rim  of  the  cylinder  is  an  absorptive  ma¬ 
terial  to  attenuate  outgoing  waves.  This  configuration  is  aimed 
to  eliminate  reflections  from  the  tissue/ambiance  interface  and 
the  cylinder/air  interface  in  order  to  enhance  the  image  resolu¬ 
tion.  An  array  of  transmitting  and  receiving  MW  antennas  are 
mounted  on  the  outer  surface  of  the  cylinder.  The  operating  fre¬ 
quency  for  the  following  examples  is  800  MHz  (except  as  oth- 
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(a)  M 


Fig.  2.  Comparison  of  the  EBA,  Bom,  CG-FFT,  and  BCGTTT  methods  and  analytical  solution  for  a  circular  cylinder  (er  =  64,  a  -  0.64  S/m)  with  the  center 
at  the  origin.  The  line  source  is  at  (0, 0.03)  m,  while  the  receivers  are  at  (0,?/)m,  y  €  {-0.01,0.01)m. 


Fig.  3.  Lossy  cylinder  with  er  —  57.2,  <r  =  1.08  S/m,  and  diameter  0.08 
m.  The  line  source  is  at  (0,0.08)m,  while  the  receivers  are  at  (0,y),y  € 
(-0.04,0.04). 

erwise  indicated),  although  both  the  apparatus  and  numerical 
models  are  designed  for  wide-hand  signals. 

A.  Forward-Modeling  Results 

In  our  2-D  forward  model,  we  simulate  an  electric  current 
source  J  =  zJz(x,  y)  in  an  inhomogeneous  conductive  medium 
with  dielectric  constant  er(x,  y),  conductivity  cr{x,  y),  and  con¬ 
stant  permeability  /xo  •  The  background  medium  is  homogeneous 
and  has  the  properties  erb  —  16,  crb  =  0.16  S/m,  and  m,  =  fi o 
to  simulate  the  normal  tissue  and  the  matching  fluid  inside  the 
cylinder. 

We  first  study  the  accuracy  of  the  EBA  method  against  the 
Bom  approximation,  CG-FFT,  BCG-FFT,  and  analytical  solu¬ 
tion  for  a  circular  tumor  with  er  —  64,  er  =  0.64  S/m  and  a 
diameter  of  0.01  m.  The  line  source  is  0.03  m  away  from  the 
cylinder  center.  Fig.  2  shows  that  the  EBA  result  has  a  much 
greater  accuracy  than  the  Bom  approximation  for  this  high  con¬ 
trast  medium.  Both  CG-FFT  and  BCG-FFT  results  have  the 
same  accuracy,  and  take  four  and  three  iterations  to  reach  a 
relative  residual  error  of  0.01%  with  the  EBA  preconditioner. 
The  CPU  time  for  the  Bom  approximation,  EBA,  CG-FFT,  and 
BCG-FFT  methods  are  1.2, 1.2,  3.5,  and  2.2  s,  respectively,  on 
a  SUN  Ultra  60  Workstation. 

For  larger  contrasts  and  higher  frequencies,  as  expected,  the 
accuracy  of  the  EBA  will  decrease,  and  full-wave  CG-  and 
BCG-FFT  methods  have  to  be  used.  In  the  second  example. 


we  simulate  a  large-contrast  circular  cylinder  with  er  —  57.2, 
cr  =  1.08  S/m,  and  a  diameter  of  0.08  m.  The  line  source  is 
0.12  m  away  from  the  center  of  the  cylinder.  We  compute  the 
electric  field  inside  the  cylinder  by  the  CG-FFT  and  analytical 
solution.  Fig.  3  is  the  comparison  of  the  numerical  result  and 
analytical  solution.  We  observe  that  both  real  and  imaginary 
parts  of  the  numerical  result  agree  well  with  the  analytical 
solution. 

Fig.  4  compares  the  CG-FFT  and  BCG-FFT  results  for  a 
tissue  with  multiple  tumors.  The  results  again  show  excellent 
agreement  between  these  two  methods.  The  BCG-FFT  method 
converges  in  only  two  iterations,  while  the  CG-FFT  method 
converges  in  eight  iterations  to  a  residual  error  of  0.01%.  The 
CPU  time  for  this  case  is  4.5  s,  and  2.2  s  for  the  CG-FFT  and 
BCG-FFT  methods,  respectively. 

B.  Inverse-Scattering  Results 

For  the  image  formation,  we  have  developed  the  two-step  in¬ 
verse  method  that  is  based  on  the  EBA  and  CSI  method.  The 
EBA  inversion  aims  to  produce  a  fast,  albeit  quantitatively  less 
accurate,  image  while  the  CSI  method  aims  to  produce  an  ex¬ 
cellent  accuracy  at  the  expense  of  more  computational  cost.  The 
two-step  inversion  result  is  used  as  the  preconditioner  for  the 
CSI  result.  In  the  following  examples,  the  measured  data  is  ob¬ 
tained  by  the  CG-FFT  method  with  a  much  finer  grid  than  that 
used  in  the  inversion.  Neither  the  EBA,  nor  the  CSI  method  re¬ 
quires  a  forward  solution  during  the  inversion  iterations,  thus, 
the  so-called  “inverse  crime”  is  not  committed  here.  The  back¬ 
ground  is  the  normal  tissue  with  properties  given  for  Fig.  2.  We 
use  60  transmitters  and  60  receivers  uniformly  distributed  on  the 
four  edges  of  a  square  of  size  16x16  cm2 ,  giving  a  total  of 3600 
measured  data  points.  The  region  to  be  inverted  is  divided  into 
29  x  29  unknown  pixels  with  a  ceil  size  0.0052  x  0.0052  m2. 

In  the  first  example,  there  are  two  anomalous  objects  within 
the  homogeneous  medium.  The  first  anomaly  has  a  dielectric 
constant  of  er  =  40  and  conductivity  cr  ~  1.08  S/m.  The  other 
object  has  er  =  57.2,  cr  —  1.08  S/m.  These  two  anomalies 
simulate  two  tumors;  the  first  one  perhaps  in  the  process  of  be¬ 
coming  malignant,  while  the  second  one  is  malignant.  There  is 
a  significant  difference  in  the  electrical  properties  of  these  two 
tumors.  The  reconstructed  results  are  shown  in  Fig.  5.  Both  the 
permittivity  and  conductivity  anomalies  are  well  resolved.  In 


LIU  et  al :  ACTIVE  MICROWAVE  IMAGING  I — 2-D  FORWARD  AND  INVERSE  SCATTERING  METHODS 


129 


Fig.  4.  CG-FFT  and  BCG-FFT  results  for  a  tissue  with  multiple  tumors.  The  line  source  is  located  at  (0,  —0.189).  (a)  Geometry  of  four  tumors  (er  =  72, 
cr  —  0.72  S/m)  in  a  background  medium  =  16,  <rb  =  0.16  S/m),  (b)  Magnitude  of  electric  field  along  {x,  0),  x  €  (-0.1, 0.1)  and  (c)  along  (0,  y),  y  G 
(—0.1, 0.1).  (d)  Convergence  curves  for  the  CG-FFT  and  BCG-FFT  methods. 

(a)  (b)  (c) 


x(cm) 


Fig.  5.  Simultaneous  inversion  of  dielectric  constant  and  conductivity  profiles  at  800  MHz.  The  ground  truth  of:  (a)  er  and  (d)  cr  profiles.  Inverted  dielectric: 
(b),  (c)  constant  and  (e),  (f)  conductivity. 

this  case,  the  preconditioned  CSI  method  takes  70  iterations.  We  then  demonstrate  the  super  resolution  [38],  [39]  of  the 
and  8  min  and  7  s  CPU  time  to  converge  to  a  relative  residual  nonlinear  inverse-scattering  CSI  algorithm.  The  same  measure- 
error  of  0.97%.  ment  geometry  as  in  Fig.  5  is  used.  We  use  an  EBA  precondi- 
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(a)  (b)  <c> 


Fig.  6.  Super  resolution  of  dielectric  constant  and  conductivity  anomalies  separated  by  only  0.22  wavelengths  at  800  MHz.  The  ground  truth  of:  (a)  er  and  (d)  o 
profiles.  The  inverted  dielectric:  (b),  (c)  constant  and  (e),  (f)  conductivity.  The  measured  data  E*zci  has  a  10%  Gaussian  noise  (SNR  =  20  dB). 


Fig.  7.  Inversion  of  dielectric  constant  and  conductivity  anomalies  separated  by  only  0.22  wavelengths  at  6  GHz  with  properties  as  in  Fig.  5.  The  ground  truth 
of:  (a)  €v  and  (d)  a  profiles.  The  inverted  dielectric:  (b),  (c)  constant  and  (e),  (f)  conductivity. 


tioned  CG-FFT  to  simulate  the  electric  field  outside  the  inver¬ 
sion  area.  Within  the  region  being  inverted,  there  are  four  anom¬ 
alies  separated  only  by  0.22  wavelengths  (2.06  cm)  at  800  MHz. 
We  add  10%  Gaussian  noise  to  the  scattered  field  data  (i.e., 
signal-to-noise  ratio  SNR  =  20  dB).  Fig.  6  shows  the  ground 


truth  and  reconstructed  profiles.  This  example  clearly  demon¬ 
strates  that  the  so-called  super  resolution  [38]  can  be  achieved 
for  the  high-contrast  medium  in  MW  breast  imaging.  This  super 
resolution  is  possible  because  the  nonlinear  multiple  scattering 
has  been  fully  accounted  for  and  because  of  the  near-field  mea- 
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surement.  For  more  in-depth  discussions  on  this  topic,  the  reader 
is  referred  to  [38]  and  [39]. 

Finally,  we  study  the  potential  resolution  at  a  higher  fre¬ 
quency  of  6  GHz.  In  this  case,  we  invert  for  two  very  small 
anomalies  only  0.28  cm  apart  and  complex  permittivity  the 
same  as  in  Fig.  5.  The  image  area  is  about  2x2  cm2  with 
29  x  29  pixels,  giving  a  comparable  computational  cost  as 
Figs.  5  and  6.  RMS  noise  of  10%  for  the  scattered  field  data 
has  been  added.  Fig.  7  shows  that  our  inverse  algorithm  can 
resolve  these  small  objects  when  operated  at  a  high  frequency. 

C.  Discussions 

This  study  represents  a  new  application  of  some  of  our  re¬ 
cent  fast  algorithms  developed  for  low-frequency  electromag¬ 
netic  subsurface  sensing  problems  [28]-[30],  [35],  [41],  [45], 
[49].  However,  for  low-frequency  induction  problems,  the  con¬ 
duction  current  is  dominant  and  only  the  conductivity  profile 
can  be  inverted.  For  MWI,  both  dielectric  constant  and  conduc¬ 
tivity  profiles  must  be  inverted  simultaneously.  This  represents 
a  much  more  difficult  problem,  and  an  increasingly  important 
application  in  biomedical  imaging. 

The  above  numerical  results  show  the  efficacy  of  our  2-D  for¬ 
ward-  and  inverse-scattering  methods.  In  spite  of  this,  there  are 
many  outstanding  issues  that  require  further  research.  Firstly, 
it  would  have  been  ideal  if  there  is  experimental  verification. 
However,  at  present,  we  are  unable  to  present  this  verification 
because  our  experimental  setup  has  3-D  waveguide  aperture 
transmitting  antennas,  which  cannot  be  adequately  modeled  by 
our  2-D  methods.  Secondly,  3-D  effects  in  breast  tumor  detec¬ 
tion  are  expected  to  be  very  strong.  Therefore,  we  are  currently 
pursuing  a  better  experimental  verification  in  the  following  two 
directions:  we  are  modifying  our  experimental  setup  so  that 
it  can  be  better  modeled  by  our  2-D  models  with  some  3-D 
compensation  techniques,  and  we  are  developing  full  3-D  for¬ 
ward-  and  inverse-scattering  methods  to  include  all  3-D  effects 
in  MWI  measurements.  Thirdly,  the  breast  tissue  is  not  homoge¬ 
neous  because  of  the  presence  of  fat,  glandular  tissue,  etc.  This 
inhomogeneity  will  definitely  degrade  the  images.  These  impor¬ 
tant  issues  will  be  addressed  in  our  future  research. 

VI.  Conclusions 

We  have  applied  fast-forward  and  inverse  methods  to  simulate 

2- D  MWI  for  breast  tumor  detection.  The  forward  methods  are 
based  on  the  EBA,  FFT,  and  CG  and  BCG  methods.  The  FFT  al¬ 
gorithm  is  used  to  expedite  the  EBA  calculation,  while  the  EBA 
is  used  to  precondition  the  CG-FFT  and  BCG  FFT  methods. 
The  combination  of  these  basic  ingredients  results  in  algorithms 
with  a  complexity  of  0{N  log2  Ar).  The  inverse  methods  are 
based  on  a  two-step  nonlinear  inversion  and  the  CSI  methods, 
which  do  not  require  forward  solutions  during  the  iterations 
of  the  inverse  solution.  Numerical  results  demonstrate  the  ef¬ 
ficiency  of  the  forward  and  inverse  methods,  and  the  high  reso¬ 
lution  of  the  image-formation  algorithm  for  the  breast-imaging 
applications.  Future  work  includes  experimental  validation  and 

3- D  algorithms. 
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ABSTRACT 

This  paper  presents  a  review  of  the  dielectric  properties  of 
normal  and  malignant  breast  tissues  for  radio  through 
microwave  frequencies,  as  well  as  a  brief  summary  of  the 
experiment  methods  and  the  mechanisms  that  explain  the 
difference  in  die  dielectric  properties  of  normal  and  malignant 
breast  tissue.  This  information  provides  a  basis  for  the 
development  of  diagnostic  techniques  for  breast  cancer  and 
also  highlights  the  areas  that  are  in  need  of  more  experiments. 

1.  INTRODUCTION 

The  contrast  in  the  dielectric  properties  between  normal  and 
malignant  tissues  is  a  basis  for  diagnostic  applications  using 
microwave  devices.  The  study  of  normal  tissues  has  been 
widely  reviewed.  This  paper,  in  addition,  collects  together 
dielectric  property  data  on  benign  and  malignant  breast  tissues 
from  a  number  of  researchers,  and  presents  them  in  graphical 
form  so  that  this  information  is  convenient  for  general 
reference.  It  should  be  emphasized  that  the  data  shown  has 
been  interpolated,  extrapolated  or  computed  from  the  graphs 
and  tables,  so  it  is  not  necessarily  precise.  This  paper  also 
reviews  the  mechanisms  behind  the  differences  in  dielectric 
properties  of  normal  and  malignant  breast  tissues. 

Most  data  are  represented  in  terms  of  conductivity  o  and 
relative  permittivity  e ,  since  cr  and  e  of  biological  materials 
are  practically  independent  of  frequency  up  to  the  microwave 
range  [1].  For  the  two  low  frequency  cases  with  no  a  and  e’ 
available,  the  data  are  represented  in  terms  of  a  parallel 
combination  of  a  conductance  G  and  a  capacitance  C.  The  two 
pairs  of  terms  are  equivalent  in  that, 

Y*  -  G  +  j©C  -  (A/d)( a+jeaeoe’)  (1) 

Where,  Y*  is  the  complex  admittance  of  the  equivalent  circuit 
of  an  idealized  parallel  plate  capacitor  filled  with  the  tissue  of 
a  and  e*.  A/d  is  the  geometry  factor.  We  assume  (e\o)  follow 
a  bi-variate  normal  distribution.  Therefore,  the  modeling  of  the 
data  can  be  represented  by  the  specific  cross  section  of  the 
distribution  function,  which  satisfies, 

{(x-mJ2/a*J+(y-iny)2/0yJ-2p(x-m]l)(y-iny)/a*cjy}/(  1  ~p2)=  1  (2) 

where,  x  is  the  relative  permittivity,  y  is  the  conductivity,  and 
m,,,  my  and  ay  are  the  marginal  mean  and  variance,  p  is  foe 
correlation  coefficient  of  e*  and  a.  If  only  the  mean  and 
variance  values  are  available  but  not  the  original  data  pairs,  we 


assume  an  independent  distribution  of  (e’,a)  and  p=0. 

All  data  are  from  human  breast  tissue,  except  for  one  case 
from  rats.  The  category  of  breast  tissues  in  the  literature  is 
ambiguous.  In  this  paper,  we  define  the  following  major 
categories  of  breast  tissue. 

•  fat 

•  normal,  includes  glands  tissue  (lobules  that  produce 
milk),  and  connective  tissue  (fibrous  tissue  that 
surrounds  foe  lobules  and  ducts) 

•  benign,  includes  fibroadenoma  and  mastitis 

•  malignant,  i,e.  breast  carcinomas 

We  first  display  the  data  in  the  low  and  high  frequency 
regions,  in  the  order  of  their  publication  date.  Then  we  discuss 
the  consistency  and  inconsistency  in  the  data  as  well  as  the 
diagnostic  value  of  foe  dielectric  properties  from  foe  data. 
Finally  the  mechanisms  are  reviewed. 

2.  REVIEW  OF  EXPERIMENT  DATA 
2.1.  List  of  data  at  loir  frequency 

1.  Fricke  etal.  ((2],  1926,  20kHz,  24°C),  measured  the  parallel 
capacitance  and  resistance  (R=l/G)  of  excised  samples  from 
55  patients,  using  a  wbeatstone  bridge.  Several  types  of  tissue 
were  studied:  fat,  gland,  mastitis,  fibroadenoma  and 
carcinoma.  Data  is  displayed  in  Fig.1.  Only  one  sample  of  fat 
is  measured,  no  variance  is  available  for  this  type. 

Figure  1:  Capacitor  versus  Resistor 
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Left:  original  data;  right:  modeling  of  foe  date 


2.  Morimoto  etal,  ([3,4], 1990,  10kHz,  37°C)  obtained  in  vivo 
measurements  of  breast  cancer,  fibroadenoma,  normal  breast 
tissue  and  Bitty  tissue  using  a  three-electrode  method.  The 
proposed  equivalent  circuit  is  composed  of  R«  parallel  with  the 
series  of  R*  and  C*.  We  transformed  it  to  parallel  R  and  C, 
using  R=RiR*/(Ri+R*),  €=0*.  With  no  original  date  and 
correlation  coefficient  available,  Fig.2  displays  the  modeling 
of  the  data,  assuming  p=0. 
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3.  Jossinet  etal.  ([5,6],  1996,  488Hz-lMHz,  21°C)  measured 
120  sanqdcs  from  64  patients,  using  impedance  probe  sensors 
connected  with  a  microcomputer  system.  In  Fig.3,  data  of  six 
types  of  tissue  is  displayed,  gland  (o),  connective  (♦),  fat  (•), 
mastopathy  (square),  fibroadenoma  (+),  carcinoma  (V).  We 
calculated  the  relative  permittivity  and  conductivity  from  the 
original  complex  impedance  data  p*  (not  the  characteristic 
impedance)  using, 

e  -  Im  (l/p>(oeo),  <* -  Rc(l/p ),  p^otjoert  (3) 


Figure  3  Permittivity  and  conductivity  versus  frequency 


2.2.  list  of  data  at  high  frequnecy 

4.  T.S.  England  etal.  ([7,8],  1949-50,  3-24GHz,  37°) 
measured  the  attenuation  a  nepers/cm,  and  phase  constant  p 
radians/cm  of  die  standing  wave  pattern  of  the  excised  human 
breast  fat  and  carcinoma  tissue  samples  in  the  wave-guide.  We 
computed  the  relative  permittivity  c  and  conductivity  0  using 
[23,Eqn4.28], 

e»-(a2-  pWnofifc  o=*2ap/©po  (4) 

The  variance  due  to  these  measurements  were  computed  and 
presented  with  the  light  marker  m  Fig.  4.  Our  extrapolations 
are  shown  with  dashed  lines. 


Figure  4  permittivity  and  conductivity  versus  frequency 


5.  W,T.  Joines  et.al.  ([9,10,11],  1980,  30MHz*2GHz,  37°C), 
obtained  in  vivo  measurements  of  SMT-2A  tumor  and 
mammary  gland  tissue  samples  from  22  rats.  The 
nondestructive  method  uses  an  open-ended  coaxial  probe  to 
produce  a  fringing  field  in  the  termination  tissue  and  a 
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Figure  5:  Permittivity  and  conductivity  versus  frequency  [9] 


6.  S.S.  Chaudhary  etal.  ([12],  1984,  3MHz-3GHz,  24°), 
measured  excised  normal  and  malignant  breast  tissues  from  1 5 
patients,  using  the  time  domain  spectroscopy  system  of  HP. 
We  use  the  total  spread  over  die  mean  value  (0.8%)  to 
compute  the  variance  of  the  data,  which  is  shown  in  Fig.  6. 


Figure  6:  Permittivity  and  conductivity  versus  frequency  [12] 

7.  AJ.  Surowiec  atel.  ([13],  1988,  0.02MHz-100MHz,  37°) 
measured  the  input  reflection  coefficient  of  28  samples  from  7 
patients,  using  a  coaxial  line  sensor  connected  to  an  HP3577 
network  analyzer.  Tissue  types  include  ductal  carcinoma, 
lobular  carcinoma,  and  surrounding  tissues.  The  measured 
dielectric  values  are  available  only  at  100kHz  and  100MHz, 
and  we  use  square  symbols  to  represent  die  mean  and 
meanistd  values  of  those  data  in  Fig.  7.  The  authors  provided 
the  parameters  e«  ,es  ^  t  o,  and  a,  by  fitting  the  data  with  the 
Ode-Cole  equations  [1], 

E*=e«+(e1-€„)/(  1  +(jfft)1'a)-j<V»ea  (5) 

where  £=l/27rr,  t  is  the  relaxation  time,  f*  is  the  relaxation 
frequency,  a  is  the  distribution  parameter  that  reflects  the 
range  of  t.  s  represents  the  low  frequencies  f  «  « 

represents  the  high  frequencies  f »  1*  When  a  *0,  Eqn.5  is 
the  same  as  Debye  equation. 


The  curves  in  Fig.  7  are  computed  from  the  fitted  Cole-Ode 
model  in  [13].  The  dark  lines  are  the  mean  values,  the  light 
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lines  are  the  mesurfcstd  values.  Extrapolations  are  shown  in 
light  color.  Comparing  the  mean  from  the  measurements  and 
die  model,  the  conductivity  values  agree  well,  but  the  relative 
permittivity  values  show  some  inconsistency. 

8.  A.M.  Campbell  eUl.  {[14],  1992, 3.2GHz,  24°Q  measured 
39  samples  of  normal  breast  fat,  IS  samples  of  benign  tumors, 
22  samples  of  glandular  connective  tissue  and  20  samples  of 
cancer  from  37  patients,  using  a  resonant  cavity  perturbation 
method.  Dielectric  properties  were  measured  using  the 
observation  of  the  changes  in  resonant  frequency.  In  Fig.8,  the 
left  plot  displays  the  original  data.  Tbe  right  plot  illustrates  the 
modeling  of  the  data  with  Eqn.2. 


Figure  8  conductivity  versus  permittivity 
Left:  original  data;  Right:  modeling  of  the  data 


9.  W.T.  Joines  ct  al.  ([15],  1994,  50MHz-900MHz,  24°Q, 
measured  admittance  of  12  normal  mammary  samples  and  12 
malignant  mammary  samples  from  12  patients,  using  a  flat- 
ended  coaxial  probe  connected  to  a  network  analyzer  HP 
8754A.  e  and  o  are  then  computed  from  the  admittance  with 
the  knowledge  of  the  geometry  factor.  In  Fig,  9,  the  mean 
values  are  presented  with  the  solid  lines,  the  standard  error  on 
the  mean  is  presented  with  the  dashed  lines.  Extrapolations  are 
represented  with  light  lines. 
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Figure  9  Permittivity  and  conductivity  versus  frequencies  [15] 


10.  P.M.  Meaney  et.aL  ([16],  2000  900MHz,37°)  obtained  the 
in  vivo  breast  microwave  imaging  of  5  patients,  all  non- 
malignant,  at  900MHz. 


Figure  10  conductivity  versus  permittivity 


We  display  the  modeling  of  the  data  of  individuals  in  Fig.  10 
to  represent  the  heterogeneity  within  and  across  patients.  The 
correlation  coefficient  is  not  available  and  assumed  to  be  0. 

23.  Data  consistency  and  inconsistency 

Low  frequencies.  At  low  frequency  ranges,  the  dielectric 
values  of  the  four  types  of  tissue  are  all  available  for  the  first 
three  cases,  as  shown  in  Figures  1,  2  and  3.  We  cannot 
compare  them  directly,  because  of  the  unknown  geometry 
factors  in  cases  1  and  2.  Yet  we  can  still  make  comparisons 
according  to  the  relative  distribution  of  the  data  for  foe  same 
tissue  types.  The  consistencies  noted  are  listed  below, 

•  The  conductivity  of  the  malignant  tissue  falls  between  the 
fat  (plus  connective  tissue)  and  the  normal  gland  tissue 
(plus  tbe  benign  fibroadenoma  and  mastitis  tissues). 

•  The  benign  and  normal  tissues  can  be  grouped  together 
relative  to  foe  malignant  tissues  on  tbe  e  -a  plane. 

•  The  relative  position  of  foe  fat  tissue  on  foe  c*-<x  plane 
compared  to  the  other  types  is  the  same. 

This  information  of  data  consistency  provides  a  basis  for 
identifying  breast  cancer,  benign  tumor  and  normal  breast 
tissue  using  the  tissue  conductivity  at  IkHz-lMHz. 

The  inconsistency  in  foe  first  three  cases  is  that  the 
malignant  tissue  has  lower  capacitance  (or  permittivity)  than 
that  of  the  normal  and  benign  tissues  in  case  2  and  the  lower 
frequency  region  of  case  3,  but  it  has  tbe  largest  capacitance 
values  in  cases  1  and  the  higher  frequency  region  of  case3. 
One  of  the  possible  reasons  for  foe  inconsistency  is  foe 
frequency  difference.  It  is  10kHz  in  case  2,  20kHz  in  case  1 
and  488Hz-lMHz  in  case  3.  Therefore,  foe  inconsistency  can 
be  related  with  a  turning  frequency  point  in  tens  of  kHz,  above 
which  foe  capacitance  value  of  cancerous  tissues  became 
larger  than  that  of  normal  and  benign  tissues.  Other  reasons  for 
the  inconsistency  are  foe  intrinsic  heterogeneity  and  the 
temperature  difference  of  tissue  samples. 

This  information  of  inconsistency  suggests  that  in  the  range 
1kHz  to  1MHz,  the  capacitance  of  breast  tissues  is  not  a  good 
quantity  to  diagnose  breast  cancer.  More  experiments  and 
analyses  on  foe  capacitance  properties  of  normal  and 
malignant  breast  tissues  are  needed  in  this  frequency  range. 

High  frequencies.  In  the  high  frequency  range,  we  can 
compare  foe  dielectric  data  directly. 


Figure  1 1  Permittivity  versus  frequency 


Proceedings  IEEE  Southeastern  2002 
459 


Figure  12  Conductivity  versus  frequency 


Fig.  1 1  and  12  illustrate  the  dielectric  data  together  in  the 
range  SOOKHz  to  20GHz.  One  more  case  of  breast  fat  [21]  is 
included.  Four  types  of  breast  tissue:  fat  {•),  normal  (o), 
benigh  (+),  malignant  (V)  are  displayed. 

For  a  clearer  view.  Fig.  13  and  14  illustrate  the  modeling  of 
the  data  for  multiple  cases  on  the  e  kt  plane  at  900MHz  and 
3.2GHz.  Again,  four  types  of  breast  tissue:  fat  (•),  normal  (o), 
benign  (+),  malignant  (V)  are  displayed. 
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Figure  13  Conductivity  versus  Permittivity,  at  900MHz 
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Figure  14  Conductivity  versus  Permittivity,  at  3.2GHz 

From  Fig.  1 1-14,  we  observe  the  data  inconsistency, 

•  The  mean  dielectric  values  of  normal  and  malignant  breast 
tissues  have  obvious  variability. 

•  The  mean  conductivity  values  of  normal  breast  tissues  in 
Joines’  rat  data,  Campbell's  data,  and  Meaney’s  data  are 
more  than  twice  that  of  other  normal  breast  tissue  cases  at 
the  corresponding  frequencies.  The  mean  permittivity 


values  of  normal  breast  tissue  in  Joines'  rat  data  and 
Campbell's  data  are  more  than  twice  that  in  England's 
data  and  Chaudhary’s  data  at  3.2GHz  frequency. 

The  possible  two  reasons  for  the  inconsistency  are  listed 
below,  which  may  help  explain  the  results  and  improve  the 
future  experiment  designs. 

•  Experiment  method 

1.  Limitations  in  experiments.  Chaudhary’s  samples 
were  collected  in  physiological  saline,  which  will 
affect  the  accuracy  of  the  data.  Campbell’s  malignant 
samples  are  frozen  and  defrosted  before  the 
measurement,  which  may  affect  the  accuracy  of  the 
data  of  this  type. 

2.  In  vivo  vs  in  vitro.  Joines*  rat  data  is  from  in  vivo 

measurements  with  uncorrected  skin  effects. 

Meaney’s  data  comes  from  reconstruction  of  in  vivo 
microwave  imaging.  Others  are  from  excised 

samples.  In  vivo  methods  seem  to  have  higher 
dielectric  values. 

3.  Sample  temperature  differences.  Lower  sample 

temperature  will  make  the  dielectric  value  a  little  bit 
lower,  when  the  frequency  is  below  2GHz. 

•  Intrinsic  heterogeneity 

1.  Normal  breast  tissues  are  composed  of  breast  fat, 
connective  tissue  and  gland  tissue,  etc.  In  the 
literature,  the  composition  of  normal  breast  tissues 
from  case  to  case  may  differ. 

2.  Different  stages  of  tumor  development  will  change 
the  tumor’s  dielectric  property  and  introduce 
variability  [2,13].  Some  samples  of  malignant  tissues 
were  actually  composed  of  small  parts  of  malignant 
cells,  infiltrating  within  a  large  part  of  normal  cells, 
which  may  decrease  the  mean  value  of  the  malignant 
tissue  samples. 

3.  Across  patients.  The  breast  tissue  samples  from 
patients  with  different  water  content  or  fet  content 
and  in  different  stage  of  menstruation,  pregnancy  or 
lactation  will  have  obvious  differences  in  dielectric 
values.  Campbell's  data  came  from  a  relatively  larger 
patient  group,  which  may  introduce  wider  variability. 

This  inconsistency  information  indicates  the  importance  of 
using  proper  sample  storage  method  before  experiments  and 
suggests  a  standardization  of  the  experiment  conditions  like 
the  sample  and  environment  temperatures,  as  well  as  the 
record  of  patients’  information  for  later  analysis. 

Although  there  are  so  many  conditions  out  of  control,  we 
still  observe  the  data  consistency  from  Fig.  1 1-14 

•  The  mean  conductivity  of  the  normal  tissue  is  less  than 
15msfcm  up  to  3.2GHz. 

•  Malignant  tissues  have  higher  mean  permittivity  and 
conductivity  values  than  those  of  normal  breast  tissues 

•  Fat  tissues  have  the  lowest  mean  permittivity  and 
conductivity  values 

This  data  consistency  information  provides  the  basis  for 
breast  cancer  diagnosis  using  the  dielectric  properties  in  the 
microwave  frequency  range. 
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24.  Discussion  of  the  diagnostic  virtues 

It  is  misleading  to  use  only  the  contrast  of  die  mean  values  to 
judge  the  diagnostic  value  of  the  dielectric  properties.  Since 
the  intrinsic  heterogeneity  in  malignant  tissue  is  large*  this  will 
decrease  the  mean  contrast  The  mean  values  from  samples 
across  patients  will  decrease  die  contrast  as  compared  to  an 
individual  patient  Therefore,  the  diagnostic  value  of  the 
dielectric  properties  seems  to  be  underestimated,  as  in  case  8. 
A  better  concept  might  be  to  use  the  contrast  of  the  maximum 
value  of  the  malignant  tissue  with  the  mean  of  the 
neighborhood  normal  tissue  samples  [2,13].  Better  criterion 
can  be  defined  using  the  probability  of  detection  and  false 
alarm,  in  which  the  random  model  of  the  point  dielectric 
values  and  the  spatially  distributed  dielectric  values  are 
incorporated. 

In  summary,  we  observed  the  diagnostic  value  of  the 
dielectric  properties  from  the  data,  as 

•  The  low  conductivity  values  of  the  normal  breast  tissue 
enable  penetration  of  microwave  frequencies  up  to  die  low 
GHz  range,  which  coincides  with  the  simulation  results  in 
[22]. 

•  At  lOOMHz-lGHz,  dielectric  properties  can  significantly 
help  classify  normal  and  malignant  tissues. 

•  At  frequency  ranges  of  lGHz-3GHz,  dielectric  properties 
can  help  classify  normal  and  malignant  tissues. 

•  At  IOkHz-lMHz,  dielectric  property  can  help  classify 
normal,  benign  and  malignant  tissues,  yet  mainly  depends 
withe  conductivity. 

2J.  Areas  in  need  of  more  experiments 

•  The  dielectric  properties  of  benign  tissues  compared  with 
that  of  the  malignant  tissues  and  normal  tissues  in  the 
frequency  range  of  100MHz-3GHz, 

•  The  spatial  distribution  of  the  dielectric  properties  of 
normal,  benign  and  malignant  breast  tissues. 

•  The  dielectric  properties  of  human  breast  cancer  in 
different  development  stages. 

3.  MECHANISM:  NORMAL  VS  MALIGNANT  TISSUE 

We  first  review  die  mechanism  of  the  dielectric  properties  of 
biological  tissues  in  general. 

The  frequency  dependence  of  dielectric  properties  of 
biological  tissues  is  related  to  the  polarization  of  molecules 
and  structural  interfaces  in  response  to  the  applied  electric 
field  [9].  Data  from  Schwan  and  Foster  on  high  water  content 
muscle  tissue  suggests  the  presence  of  three  dispersion 
regions:  alpha,  beta,  and  gamma,  with  the  relaxation 
frequencies  to  be  KHz,  hundreds  of  kHz,  and  GHz  [17].  The 
delta  dispersion,  located  in  half  way  between  beta  and  gamma 
regions,  has  also  been  identified  [1,18,19]. 

For  engineering  applications,  the  alpha  dispersion  has  little 
significance  [1].  Beta  dispersion  occurs  at  radio  frequencies, 
and  arises  principally  from  the  charging  of  cellular 
membranes,  with  smaller  contributions  from  the  protein 
constituents  and  ionic  diffusion  along  surfaces  in  the  tissue 
[1,18,19], 


Tissues  typically  exhibit  a  small  dispersion  between  0.1  and 
3  GHz,  which  have  been  termed  the  delta  dispersion  [1,19]  or 
“UHF  relaxation”  [18].  A  combination  of  mechanisms  are 
suggested  for  this  region:  bipolar  relaxation  of  the  water  of 
hydration  “bound”  to  proteins,  a  Maxwell-Wager  effect  due  to 
ions  in  the  cytopolsm  collection  against  relative  nonconductive 
protein  surface  and  rotation  of  polar  side-chains  on  the  protein 
surface  [1,9,18],  The  relaxation  frequency  is  dominant  mostly 
by  bound  water  (ft  of  100-1000  MHz  [1],  Protein  molecules 
(ft  of 40-300  MHz)  and  free  water  (/r=25GHz)1 
The  gamma  dispersion  occurs  with  a  center  frequency  near 
25GHz  at  body  temperature,  due  to  the  dipolar  relaxation  of 
the  free  and  bound  water  and  ionic  conductivity.  Campbell  and 
Land  [14]  attribute  higher  than  expected  conductivity  at  3.2 
GHz  to  the  “tail  end”  of  fi-dispersion  effects. 

In  recent  studies,  a  variety  of  factors  have  been  explored, 
which  lead  to  pronounced  difference  in  dielectric  properties  in 
normal  compared  with  malignant  tissues,  as  listed  below: 

Necrosis.  Inflammation  and  necrosis  are  commonly  found  in 
malignant  breast  tissues.  Presence  of  necrosis  leads  to 
breakdown  of  cell  membranes  and  thus  a  larger  fraction  of  the 
tissue  that  can  carry  current  at  low  frequencies  [20],  which 
decreases  die  capacitance  of  the  tumor  [2]. 

Charging  of  the  cell  membrane.  In  breast  carcinoma,  there  is 
a  progressive  replacement  of  fat  lobules  with  fibroblastic 
proliferation  and  epithelial  cells.  Which  also  accompanied  by  a 
variety  of  alterations  at  the  transformed  cell  surface  [12(22)]. 
Cancer  cells  have  reduced  membrane  potentials  and  tend  to 
have  altered  ability  to  absorb  positive  ions  [19(77)],  they  have 
a  higher  negative  surface  charge  on  their  membranes 
[9(4,21),  19(78)].  According  to  Joines  et.alM  conductivity  of 
the  malignant  tissues  is  increased  with  this  mobile  charge 
being  displaced  and  rotated  by  the  microwave  field  [9]. 

Relaxation  tunes.  The  relaxation  times  in  malignant  tissues 
are  larger  than  those  in  normal  tissue,  indicating  that  a 
significant  increase  in  the  motional  freedom  of  water  has 
occurred  [19(75)].  Surowiec  etal.[13]  reported  that  cancerous 
breast  tissues  have  average  dielectric  relaxation  times  between 
0.6us  and  1.4  us  and  the  surrounding  normal  tissues  had 
shorter  relaxation  times  of  0.3us. 

Sodium  concentration  and  water  content.  The  sodium 
concentration  in  tumor  cells  is  higher  than  in  normal  cells 
[19(76)].  The  excessive  sodium  concentrations  not  only  affect 
the  cell  membrane  potentials  [10,  (9-1 1)],  but  causes 
malignant  tissue  to  retain  more  fluid.  According  to  Joines  et.al, 
the  excess  sodium  fluid  alone  would  yield  greater  conductivity 
and  permittivity  values  in  malignant  tissue  than  in  normal 
tissue.  In  addition,  the  fluid  is  retained  in  the  form  of  bound 
water,  which  has  larger  values  of  a  and  e  than  free  water  [10]. 

Malignant  tissues  have  significantly  higher  water  contents 
than  normal  tissues  [18,19  (76)].  The  data  from  Campbell  and 
Land  [14]  illustrates  the  dielectric  properties  related  with  tiie 
water  content  at  3. 2 GHz  of  the  breast  tissues,  as  shown  in  Fig. 
15  and  16,  The  relationship  between  relative  permittivity  and 
water  content  is  strikingly  similar  to  the  relationship  between 
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conductivity  and  water  content  This  leads  to  the  conclusion 
that  the  same  mechanism  is  responsible  for  the  change  in  both 
dielectric  properties. 

Malignant  breast  tissue  has  a  higher  ratio  of  water  content 
compared  with  that  of  the  normal  tissue,  which  coincides  its 
higher  values  of  permittivity  and  conductivity  than  normal 
breast  tissue  at  the  same  microwave  frequency.  However,  in 
this  data  taken  at  3.2  GHz,  there  is  not  a  marked  difference  in 
the  water  content  of  benign  breast  tissue  and  malignant  tumor. 


20  30  40  50  60  70  80 

Water  Content  (%  by  weight) 

Figure  1 5  Relative  permittivity  of  human  breast  tissue  vs 
water  content 
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Water  Content  (%  by  weight) 

Figure  16  Conductivity  of  human  breast  tissue  vs.  water 
content 

4.  CONCLUSION 

This  paper  presents  an  initial  review  and  consolidation  of  the 
dielectric  properties  of  normal,  benign  and  malignant  tissues  in 
the  range  of  10kHz-20GHz.  A  brief  explanation  of  the 
experiment  methods  is  presented  as  well  as  the  mechanisms 
that  explain  the  difference  in  the  dielectric  properties  of 
normal  and  malignant  tissues.  The  consistency  and 
inconsistency  of  the  data  are  discussed  as  well  as  suggestions 
for  the  possible  inconsistency.  It  is  observed  that  the  dielectric 
properties  of  breast  tissue,  even  though  containing  uncertainty, 
have  good  diagnostic  value  in  the  range  of  100MHz-3GHz. 
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ABSTRACT 

In  this  paper  the  Markov  Random  Field  is  used  to  model  the 
breast  permittivity  cross  section  as  a  propagating  medium, 
and  incorporate  it  into  the  forward  Electromagnetic  (EM) 
propagation  to  predict  the  random  field  of  the  EM  measure¬ 
ments  at  a  received  array  of  sensors.  Given  these  EM  field 
measurements,  Bayesian  approaches  are  then  developed  to 
compute  the  likelihood  ratio  for  tumor  detection  and  the  a 
posteriori  probability  display  of  tumor  localization.  Quan¬ 
titative  performance  evaluations  using  simulations  demon¬ 
strate  the  advantage  of  using  the  Bayesian  approach  to  di¬ 
rectly  process  the  measurement  data  as  compared  to  using 
the  Bayesian  or  threshold  approaches  to  detect  and  localize 
the  tumor  based  on  the  reconstructed  permittivity  image. 

1.  INTRODUCTION 

Microwave  imaging  is  a  promising  sew  modality  for  breast 
cancer  diagnosis,  partly  because  it  is  non-invasivc  and  the 
permittivity  contrast  between  normal  and  malignant  breast 
tissues  is  Ugh  [11.  In  addition,  the  attenuation  of  EM  propa¬ 
gation  in  normal  breast  tissues  is  low  so  that  it  can  penetrate 
into  the  depth  of  die  tissue  1 1].  Most  of  the  research  in  this 
field  has  focused  on  the  study  of  the  dielectric  properties 
[l,(l-22)],  the  design  of  the  microwave  imaging  prototypes 
[21,  and  the  improvement  of  the  EM  forward  and  inverse 
algorithms  [3][4][5].  However,  none  of  this  research  has 
incorporated  signal  detection  theory  directly  into  the  mi¬ 
crowave  imaging  at  the  measurement  teveL  Markov  Ran¬ 
dom  Fields  (MRF  Hammersley  and  Clifford  [6])  and  detec¬ 
tion  theory  have  been  applied  in  mammagraphy  for  diagno¬ 
sis,  such  as  [7][8][9J.  However,  [7]  and  [8]  only  assumed 
a  simple  deterministic  disk  object  model.  [9]  did  not  uti¬ 
lize  the  a  priori  knowledge  of  the  projections.  This  paper 
presents  Bayesian  algorithms  for  Scattered  Electromagnetic 
fields  through  an  Uncertain  Permittivity  Image  ( BPJSEUPl ), 
which  incorporates  the  knowledge  of  the  a  prior  permittiv¬ 
ity  image  modeled  by  the  MRF,  the  measurement  noise,  as 
well  as  the  physical  model  of  the  forward  scattered  elec¬ 
tric  field.  The  Bayesian  algorithms  for  the  Uncertain  Per¬ 
mittivity  Image  (BPJJPl)  and  the  Threshold  Image  Proces¬ 


sors  (TIP)  are  also  presented  for  comparisons.  It  should  be 
noted  that  the  forward  EM  scattering  field  is  computed  us¬ 
ing  the  Extended  Born  Approximation  (EBA)  accelerated 
CGFFT  method,  which  has  been  proposed  by  Zhang  and 
Liu  in  [41(10].  The  reconstructed  permittivity  image  is  ob¬ 
tained  using  the  EBA  as  the  initial  solution  followed  by  the 
Contrast  Source  Inversion  (CSI  [11])  method,  which  has 
been  proposed  by  Zhang  and  Liu  [4]. 

2.  DIAGNOSIS  MODEL 

The  binary  hypotheses  considered  are: 

Ho:  No  tumor  present 

Hi:  Thmortsize  LxL)  present,  located  at  an  unknown  posi¬ 
tion  S  on  the  2D  lattice  H  of  the  permittivity  cross  section. 

The  decision  as  to  whether  the  tumor  is  present  or  not, 
I.e.  whether  Ho  or  Hi  is  true,  is  the  detection  problem. 
Where  the  tumor  is  located  if  Hi  is  true,  is  considered  to  be 
a  localization  problem. 

The  data  for  the  BP-SEUPI  arc  the  microwave  measure¬ 
ments  r,  shown  in  the  middle  of  the  Fig.  1-  r  is  composed  of 
concatenating  pieces.  Each  piece  is  a  complex  Kx  1  vector 
representing  the  narrow-hand  frequency  component  of  the 
scattered  electric  field  sampled  by  K  sensors,  from  a  sin¬ 
gle  transmitter,  which  is  one  of  the  K  sensors.  It  assumes 
r-$  +  n,sisthe  signal,  and  n  is  the  additive  noise  at 
the  sensors,  modeled  by  a  multivariate  complex  Gaussian 
distribution,  with  zero  mean  and  o^Ik?  covariance  matrix. 
The  signal  to  noise  ratio  (SNR)  is  given  by  10lo£lQ(kf~) 
(dB).  The  propagating  medium  of  the  EM  field  is  the  un¬ 
certain  permittivity  image  modeled  by  the  MRF,  a  sample 
of  which  is  shown  in  the  left  of  Fig.  I. 

The  data  for  BPJJPl  and  TIP  is  the  reconstructed  per¬ 
mittivity  image  er  =  as  shown  in  the  right  of  Fig. 

1  where  T-1  represents  the  reconstruction  procedure. 

The  uncertainties  in  this  problem  are  the  tumor  position, 
the  noise  at  the  sensors,  as  well  as  the  spatial  distribution  of 
the  breast  permittivity  which  reflects  statistically  the  tissue 
variance  in  individuals  and  the  background  structure  vari¬ 
ance  across  patients. 
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Fig.  1.  Numbers  1-24  represent  the  sensors.  Left:  original 
tissue  permittivity  image;  Middle:  real  and  imaginary  part 
of  a  sample  measurement  data;  Right:  reconstructed  tissue 
permittivity  image. 

2.1.  Model  of  Breast  Permittivity  Image 

We  utilize  the  MRF  to  model  the  breast  permittivity  im¬ 
age,  because  it  implements  the  idea  that  spatially  nearby  tis¬ 
sue  permittivities  are  similar.  The  Gaussian  MRF  (GMRF, 
Chellappa  [6])  is  selected  because  it  provides  a  simple  form 
to  be  incorporated  to  the  Bay  sian  approaches.  It  is  important 
to  realize  that  this  model  is  statistical,  and  is  not  intended  to 
be  a  model  of  the  detailed  deterministic  structure  of  permit¬ 
tivity  of  a  particular  individual.  The  GMRF  is  defined  as 

«*p  flwMw»(*w-w*Wa>  ^ 

1  * 


(T.G.Birdsall).  The  optimal  localization  processor  calcu¬ 
lates  the  a  posteriori  probability  of  the  tumor  position  given 
the  data  vector.  We  derive  the  likelihood  ratios(A)  for  the 
Bayesian  detector  and  threshold  detector,  as  well  as  the  a 
posteriori  probability  image  for  the  Bayesian  localization 
processors.  We  assume  that  the  unknown  tumor  position 
has  a  uniform  distribution  a  priori  on  the  2D  lattice  fl. 


2.2.1.  Bayesian  Processor  for  Scattered  EM  field  propa¬ 
gated  through  theUncertain  Permittivity  Image  (BP SEUPI) 

detector  A(r)  oc 

Ssen  exp  (-  )p(er  ]  gl,S)&r 

f(  exp(-<r-T<^-'rC^)))p(£r  |  go )der 

(2) 

localization  processor 


P(S|r)  oc/eexp(-4 


l,S)<fcr. 
(3) 

Both  the  BP,SEUPT  detector  and  the  localization  proces¬ 


sor  require  a  high  dimensional  integration,  which  is  a  tough 
problem.  This  paper  tries  to  overcome  this  difficulty  by  us¬ 
ing  a  multivariate  complex  Gaussian  distribution  to  approx¬ 
imate  the  random  field  of  the  measurements  data,  which  has 
been  found  to  have  good  performance  in  the  simulations. 


where  and  are  the  mean  and  variance  at  (ij).  £(#),(«) 

are  the  interaction  coefficients.  The  P's  provide  two  kinds 
of  information:  one  is  the  neighborhood  of  (i  j),  i.e. 
which  implies  that  the  permittivity  at  (i  j)  only  depends  on 
the  permittivities  in  it’s  neighborhood  Another  is  the  influ¬ 
ence  of  the  neighbors  on  the  point  (ij),  as  represented  by 
the  sign  and  value  of  the  interaction  coefficients. 

We  assume  the  permittivity  image  size  is  9.2cmx9.2cm, 
529  pixels.  The  tumor  size  is  25  pixels.  The  mean  values  of 
tumor  and  background  are  pt  =  40  and  =  30,  whose  val¬ 
ues  are  similiar  to  normal  and  malignant  breast  permittivity 
values  from  experiments  [1].  The  interaction  coefficients 
of  tumor  and  background  are  fit  with  3  pixel  correlation 
length  and  ft,  with  20  correlation  length.  An  algorithm  pro¬ 
posed  by  Rue  [12]  is  used  to  fit  the  interaction  coefficients 
to  the  Gaussian  field.  Applying  Hammersley-Clifford  theo¬ 
rem  [6],  given  Eqn.l,  we  derive  the  mean  vector  p  and  ps, 
covariance  matrix  Q  and  Qs  of  the  joint  pdf  of  the  permit¬ 
tivity  image  under  the  Hq  and  Hi  conditions  respectively. 
Subscript  S  denotes  the  tumor  position. 

2.2.  Detection  and  Localization  Approaches 

According  to  signal  detection  theory,  the  optimal  detector  is 
the  likelihood  ratio  of  the  data  vector  followed  by  a  thresh¬ 
old  whose  value  is  determined  by  the  optimum  criterion 


2.2.2,  Bayesian  Processor  for  Uncertain  Permittivity  Im¬ 
age  (BP.UPI) 

detector  A(er)  oc 

Es€n  I  Qs  I 

(4) 

localization  processor 

p(S|e,)  oc|  Qs  I  e*P  (-  ~  ~  )  <5> 

2.2.3.  Threshold  Image  PwcessorfTIP) 

detector  A  =  max  tr  (?) 

localization  processor  S  =*  «t(5)  (8) 

3.  SIMULATION  RESULTS 
3.1.  An  example 

Fig.  2a)  shows  an  example  of  a  stochastic  background  per¬ 
mittivity  image  of  the  tissue,  along  with  a  simulated  tumor, 
modeled  by  the  GMRF.  Fig.  2b>d)  are  the  reconstructed 
permittivity  images  from  the  perfect  measurement  data  as 
well  as  from  60dB  and  50dB  noisy  measurement  data.  The 
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signal  detection  approach  using  the  permittivity  data  com¬ 
putes  the  a  posteriori  probability  of  the  tumor  location  given 
either  the  original  permittivity  image,  or  the  reconstructed 
image  as  data.  Fig.  2e  gives  an  upper  bound  on  tumor  lo¬ 
calization  by  plotting  the  a  posteriori  probability  of  tumor 
location  using  the  tissue  data  of  Fig.  2a.  Fig.  2f-h  shows  the 
a  posteriori  plots  of  tumor  location  based  on  post  process¬ 
ing  the  reconstructed  tissue  data  shown  in  Fig.  2b-d.  Rg. 
2j-l  shows  the  a  posteriori  plot  of  tumor  location  based  on 
the  same  measurements  used  to  get  the  reconstructions  in 
Rg.  2b-d. 

Plots  2g-h  show  that  at  the  60dB  and  50 dB  SNR  con¬ 
dition,  the  BP_UPI  using  the  reconstructed  data  misses  the 
correct  location  of  the  tumor,  and  Rg.  2k-i  shows  that  the 
BP_SEUPI  using  the  measurement  data  gets  the  correct  tu¬ 
mor  localization  with  high  probability.  This  is  a  specific  ex¬ 
ample  where  the  BP-SEUPI  works  better.  In  the  following 
sections,  it  is  demonstrated  statistically  that  the  processors 
using  the  measurement  data  have  better  performance. 


Fig.  2.  a)  Original  tissue  permittivity  image;  b)-d)  Recon¬ 
structed  permittivity  from  the  measurement  data  b)  with¬ 
out  additive  noise  c)d)  with  additive  noise,  c)  60dB  SNR  d) 
50dB  SNR;  e)-h)  the  a  posteriori  probability  of  the  tumor 
position  given  the  permittivity  image  data  -  p(S  |  er),  data 
er  comes  from  a)-d);  j>l)  the  a  posteriori  probability  of  the 
tumor  position  given  the  measurement  data  -  p(S  |  r) 


3.2.  Detection  performance 

Fig.  3a)  illustrates  the  detection  performance  comparisons 
assuming  no  additive  noise.  In  3a)  tire  ROC  of  the  thresh- 
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Fig.  3.  The  detection  performance  of  the  BP-UFI  using 
original  tissue  data,  BP-UFI  using  reconstructed  permittiv¬ 
ity  image  data,  IIP  using  reconstructed  permittivity  image 
data  and  BP-SEUPI  using  measurement  data.  a)no  additive 
noise  presents  at  sensors;  b)c)  [Inf,  60,50]  dB  noise  at  sen¬ 
sors  b)  BPJLJPI  &  TIP  c)  BP-UPI  &  BP-SEUPI 


old  detector  provides  a  performance  lower  bound.  The  TIP 
mimics  the  way  a  routine  visual  examination  of  the  image 
might  be  done.  Although  it  is  not  sophisticated,  it’s  ROC  re¬ 
flects  the  problem  of  high  positive  predictive  VaIue(PPY)  of 
conventional  mammography.  Using  the  reconstructed  per¬ 
mittivity  image  data,  the  BP-UPI  detector  is  much  better 
than  the  TIP  detector,  especially  when  the  probability  of 
false  alarm  is  low,  because  the  BP-UPI  utilizes  the  a  pri¬ 
ori  knowledge  of  the  tissue  background  across  patients  and 
the  a  prior  knowledge  of  the  different  characteristics  of  the 
normal  and  malignant  tissues  to  improve  the  detection  per¬ 
formance. 

Rg.  3a)  also  shows  that  the  BP-SEUPI  detector  using 
the  measurement  data  is  better  than  BP-UPI  detector  using 
the  reconstructed  permittivity  data,  yet  worse  than  BP-UPI 
detector  using  the  original  permittivity  data.  In  reality,  we 
do  not  have  access  to  the  original  tissue  permittivity  im¬ 
age  directly  but  to  the  EM  measurements.  However,  this 
provides  an  upper  bound  for  performance  evaluations.  The 
BP-UPI  detector  using  the  original  permittivity  data  is  better 
than  the  BP-SEUPI  because  the  forward  EM  field  maps  the 
variables  from  the  original  permittivity  domain  to  the  mea¬ 
surement  domain,  which  shrinks  the  random  variable  space 
and  decreases  the  detectablity. 

Rg.  3b)c)  compares  the  performance  of  the  detectors 
when  noise  is  present.  It  indicates  that  sensor  noise  de¬ 
grades  all  die  detection  performances.  It  also  demonstrates 
that  for  three  SNR  conditions,  both  Bayesian  detectors  are 
better  than  the  threshold  detector  and  the  BP-SEUPI  is  bet¬ 
ter  than  the  BP-UPI  using  the  reconstructed  permittivity  data. 
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3*3.  Localization  performance 

Fig.  4  shows  the  localization  performance  of  the  BP.SEUPI 
using  the  measurement  data,  and  the  BP.UPI  and  TIP  using 
the  reconstructed  tissue  permittivity  data.  The  localization 


Fig.  4.  Localization  performances:  solid  line,  BPJSEUPI 
using  the  measurement  data;  dotted  line,  BP.UPI  using  the 
reconstructed  tissue  permittivity  data;  dashed  line,  TIP  us¬ 
ing  the  reconstructed  tissue  permittivity  data 

performance  is  shown  using  the  probability  of  correct  lo¬ 
calization  (PCL)  curves.  PCL  is  obtained  by  computing  the 
ratio  of  the  number  of  correct  localizations  over  the  total 
trials.  The  localization  is  correct  if  the  located  tumor  po¬ 
sition  is  within  die  test  window.  If  it  is  requited  that  the 
located  position  is  exactly  the  same  as  the  real  position  to 

be  true,  the  window  diameter  is  zero.  For  the  other  values 
of  the  window  diameter  1,2,2.8284,  the  window  sizes  are 
9,13,25  pixels.  Fig.  4  shows  that  at  50dB  or  6GdB  SNR 
condition,  both  the  BP.UPI  and  the  BPJSEUPI  have  better 
performance  than  0.9  PCL.  They  arc  much  better  than  the 
threshold  approach  for  all  the  tested  SNR  conditions.  The 
BPJSEUPI  localization  using  the  measurement  data  is  the 
best  of  the  three. 


4w  CONCLUSION 

The  results  of  the  detection  and  localization  performances 
show  that  the  decision-aided  Bayesian  microwave  imaging 
approach  proposed  in  this  paper,  BP-SEUPI,  has  the  poten¬ 
tial  of  providing  additional  and  useful  information  for  radi¬ 
ologists.  It  is  an  algorithm  whose  goal  is  at  the  heart  of  the 
ultimate  decisions;  i.e.  is  a  tumor  present,  and  if  so  where. 
It  benefits  from  incorporating  the  a  priori  knowledge,  al¬ 
though  uncertain,  of  normal  and  malignant  breast  permit 
tivity.  It  also  gains  a  decision  performance  advantage  by 
processing  the  measurements  directly  rather  than  as  a  post 
processor  to  a  reconstructed  image. 
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Incorporating  statistical  models  of  tissue 
permittivity  into  an  optimal  signal  detection  theory 
framework  for  the  detection  of  breast  cancer  at 

microwave  frequencies 

Liewei  Sha,  Loren  W.  Nolte,  Zhong  Qing  Zhang  and  Qing  Huo  Liu 


Abstract —  This  paper  addresses  directly  the  decision- 
theoretic  task  of  detection  and  localization  of  breast  tumor 
microwave  measurements  and  optimal  signal  detec¬ 
tion  theory.  Microwave  energy  has  the  advantages  that  at 
low  power  levels  there  are  no  radiation  dangers,  no  con¬ 
trast  agents,  and  the  examinations  are  comfortable.  An  op¬ 
timal  signal  detection  theory  framework  is  used  to  improve 
the  probability  of  tumor  detection  for  a  fixed  probability 
of  false  alarm.  2D  Markov  Random  field  statistical  models 
that  characterize  the  spatial  properties  of  both  benign  and 
malignant  breast  tissue  permittivity  are  incorporated  into 
the  forward  electromagnetic  (EM)  propagation  to  predict 
the  random  field  of  the  EM  measurements  at  a  received  ar¬ 
ray  of  sensors.  This  Information  in  turn  is  incorporated  into 
a  likelihood  ratio  algorithm  for  tumor  detection.  Examples 
of  Receiver  operating  characteristics  (ROC)  and  probabil¬ 
ity  of  correct  localization  (PCL)  curves  are  presented  as 
a  function  of  local  uncertain  tissue  permittivity  character¬ 
istics,  and  tumor  contrast,  size,  and  shape,  and  compared 
with  algorithms  that  optimally  post  process  a  reconstructed 
image.  Computationally  simpler  tumor  detection  and  local¬ 
ization  algorithms,  simulations  of  their  convergence,  and  ap¬ 
proximations  to  the  ROC  using  the  detectability  index  are 
presented.  Simulations  also  indicate  the  effect  of  the  sensor 
array  configuration  on  tumor  detection  performance. 

Keywords —  Bayesian,  detection,  localization,  microwave 
imaging,  breast  tumor,  permittivity 

I.  Introduction 

Breast  cancer  is  a  significant  public  health  problem  for 
women  in  the  world.  The  major  traditional  modality  for 
breast  cancer  diagnosis  is  X-ray  mammography.  It  is  rel¬ 
atively  cheap  and  fast,  but  it  exposes  the  body  to  ionized 
radiation.  In  order  to  operate  within  safe  limits,  the  con¬ 
trast  of  the  images  tends  to  be  low.  Microwave  imaging 
is  a  promising  new  modality  for  breast  cancer  diagnosis, 
because  it  is  non-invasive  and  the  normal  and  malignant 
tissues  of  the  breast  have  high  contrast  in  dielectric  prop¬ 
erties  in  certain  ranges  of  the  electromagnetic  (EM)  fre¬ 
quencies^],  [2],  [3],  [4],  [5],  [6].  In  addition,  the  attenua¬ 
tion  of  EM  propagation  in  normal  breast  tissues  is  low  so 
that  it  can  also  penetrate  into  the  depth  of  the  tissue[7], 
[8],  [9],  [10].  At  this  time,  breast  microwave  imaging  is  not 
intended  to  replace  X-rays,  but  to  provide  additional  in¬ 
formation  to  the  radiologist  to  improve  the  performance  of 
breast  cancer  diagnosis. 

The  Authors  are  with  the  Department  of  Electrical  and  Computer 
Engineering,  Duke  University,  Durham,  NC  27708-0291  USA.  E-mail: 
lwn@ee.duke.edu 

Submitted  to  the  IEEE  Transactions  on  Medical  Imaging 


Most  of  the  research  in  this  field  has  focused  on  the  ex¬ 
perimental  and  theoretical  study  of  the  dielectric  properties 
of  breast  tissue  at  microwave  frequencies^- 10],  the  design 
of  microwave  imaging  prototypes  [11]  and  the  improvement 
of  2D  and  3D  EM  reconstruction  algorithms  [12], [13],  and 
[14].  However,  none  of  this  research  has  incorporated  signal 
detection  theory  directly  into  microwave  imaging.  Bayesian 
theory  has  been  applied  in  mammography  for  diagnosis, 
such  as  [15]  [16]  [17].  However,  [15]  and  [16]  did  not  incor¬ 
porate  a  random  model  of  the  tissue  being  imaged,  but 
assumed  a  simple  deterministic  disk  object  model.  In  [17] 
the  reconstructed  tissue  image  was  post-processed,  but  the 
propagation  model  was  not  fully  utilized.  In  [18],  initial 
performance  analysis  was  presented  using  a  Bayesian  mi¬ 
crowave  imaging  approach.  This  paper  includes  the  com¬ 
parison  of  two  Bayesian  signal  detection  theory  approaches 
for  tumor  detection  and  localization.  In  one  approach  the 
spatial  uncertainty  of  the  tissue  permittivity  is  included  in 
the  forward  EM  propagation,  and  optimal  signal  detection 
theory  is  applied  to  the  raw  received  sensor  information,  i.e. 
without  prior  reconstruction  of  an  image.  In  the  other  ap¬ 
proach,  signal  detection  theory  is  applied  as  a  postproces¬ 
sor  to  the  reconstructed  image. 

In  this  paper  the  physics-based  statistical  decision  the¬ 
ory  model  is  introduced  first.  Second,  the  Markov  Random 
Field  (MRF)  [19]  is  presented  which  is  used  to  represent 
the  uncertainty  in  the  permittivity  of  the  breast  tissue. 
Optimal  physics-based  signal  detection  theory  is  used  to 
develop  an  image  processing  algorithm  for  both  detection 
and  tumor  localization  that  incorporates  the  uncertainty 
of  EM  propagation  through  a  tissue  medium  with  uncer¬ 
tain  permittivity.  Optimal  signal  detection  theory  is  also 
used  to  develop  an  image  processing  algorithm  that  incor¬ 
porates  uncertainty  as  a  post  processor  to  a  reconstructed 
image.  Algorithms  are  also  developed  for  obtaining  upper 
and  lower  bounds  of  detection  performance,  some  of  which 
greatly  improve  the  speed  with  which  performance  bounds 
can  be  computed.  Finally,  simulation  results  are  presented 
for:  1)  the  detection  and  localization  performance  as  a 
function  of  tumor  contrast,  size,  local  uncertain  tissue  per¬ 
mittivity  characteristics,  and  shape.  2)  the  convergence  of 
the  model  covariance  matrix  estimation  as  support  for  the 
number  of  Monte  Carlo  trials  for  the  Bayesian  approaches 
3) comparison  of  the  optimal  physics-based  detection  and 
localization  approaches  that  incorporate  uncertainty  in  the 
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EM  propagation  with  reconstructed  image  post  processing 
and  4)  the  effect  of  the  sensor  array  configuration  on  the 
performances.  The  appendix  includes  a  proof  of  the  in¬ 
equality  of  the  signal  detection  theory  deflection  coefficients 
from  original,  measurement  and  reconstruction  data. 


II.  Physics-based  statistical  decision  theory 

MODEL 

In  this  paper  the  problem  of  malignant  tissue  detection 
is  placed  within  the  framework  of  physics-based  driven  sig¬ 
nal  detection  theory.  The  decisions  are  first,  the  detection 
problem,  is  a  tumor  present  or  not  somewhere  in  the  tissue? 
Secondly,  the  localization  problem,  if  present  somewhere, 
where  is  its  location?  We  are  particularly  interested  in  im¬ 
proving  the  early  detection  and  localization  performance 
for  the  more  difficult,  and  less  obvious,  situations  due  to 
uncertainties  in  tissue  characteristics.  The  detection  prob¬ 
lem  can  be  expressed  by  the  binary  hypotheses 
Hi:  Malignant  tissue  present,  regardless  of  its  particular 
location,  shape,  and  microwave  permittivity  characteris¬ 
tics, 
or 

Hq:  No  malignant  tissue  present 

The  basic  objective  is  to  make  optimal  decisions  about 
whether  Hi  or  H0  is  true  and  if  Hi  is  true,  what  is  the 
best  estimate  of  the  size,  location,  and  permittivity  char¬ 
acteristics  of  the  tissue. 

The  GMRJF  model  presented  in  section  Illis  used  to  char¬ 
acterize  the  uncertainty  in  the  permittivity  of  benign  and 
malignant  tissue.  An  example  realization  of  a  2D  cross 
section  of  tissue  permittivity  using  this  model,  and  a  par¬ 
ticular  sensor  array  configuration,  is  illustrated  in  Figure 
la.  In  this  case,  24  sensors  are  arranged  in  a  rectangle  sur¬ 
rounding  the  tissue  It  is  important  to  realize  that  our  char¬ 
acterization  is  statistical  so  that  the  model  does  not  convey 
detailed  deterministic  anatomical  features.  The  conductiv¬ 
ity  data  is  assumed  to  be  a  known  constant  for  simplicity. 
The  numbers  1-24  show  the  positions  of  sensors  of  one  of 
the  sensor  array  configurations  used  in  this  paper.  All  24 
sensors  can  both  transmit  and  receive  microwave  signals, 
as  in  a  tomographic  situation. 

The  microwave  measurements  at  the  sensors  are  repre¬ 
sented  by  a  K2  x  1  complex  row  vector  r  =  [ri,r2...rit], 
where  r»  is  a  K  x  1  complex  row  vector,  representing  the 
data  relevant  to  the  ith  transmitter.  K  is  the  number  of 
sensors.  An  illustration  of  one  realization  of  the  measure¬ 
ment  vector  r,  real  and  imaginary  components,  is  shown  in 
figure  lb,  with  K2  —  576.  These  measurements  are  result 
of  the  scattered  electromagnetic  field  that  has  propagated 
through  the  uncertain  permittivity  media.  This  is  an  ex¬ 
ample  of  data  used  by  the  Bayesian  approach  which  incor¬ 
porates  directly  the  uncertainty  of  the  raw  measurements. 
More  specifically,  under  the  H\  hypothesis,  the  data  vector 


a) 


Fig.  1.  a)  original  tissue  permittivity  image,  b)  measurements  of 
scattered  electric  field  (mV/m),  c)  reconstructed  tissue  permittivity 
image.  Numbers  1-24  represent  the  array  configuration. 

r  can  be  expressed  as: 
r  —  s  +  n 

Si  =  E^ip^qr),  R  =  1...K,  T  =  1...K 

SNR(dB)  =  1^  =  ^ 

(1) 

where  R  represents  the  sensor  as  a  receiver  and  T  represents 
the  sensor  as  a  transmitter.  By  using  each  sensor  as  the 
transmitter  in  turn,  and  all  sensors  as  receivers,  we  obtain 
the  data  vector  r.  s  is  the  pure  signal,  a  K2  x  1  complex 
vector,  n  represents  the  additive  noise  at  the  sensors,  n  ~ 
Af(0, 2 o%Ik*)-  AT  notifies  a  multivariate  complex  Gaussian 
distribution  ([20]).  I#?  is  an  Identity  covariance  matrix 
with  size  K 2  x  K 2. 

The  other  data  representation  that  is  used  is  that  of  the 
reconstructed  permittivity  image  er  =  T_1(r),  an  example 
of  which  is  shown  in  figure  1c,  where  T~l  represents  the 
reconstruction  procedure. 

Many  different  types  of  uncertainties  can  be  incorporated 
into  the  information  processing.  In  this  paper,  the  princi¬ 
pal  uncertainties  considered  are  the  the  spatial  uncertainty 
of  both  benign  and  malignant  breast  tissue  permittivity, 
tumor  position,  its  size  and  shape,  and  the  noise  at  the 
sensors. 

A.  Measurements 

The  pure  measurement  data  is  obtained  by  computing 
the  scattered  electric  field  through  the  uncertain  tissue  me¬ 
dia,  using  the  Extended  Boro  Approximation  (EBA)  im¬ 
plemented  with  the  improved  CGFFT  method,  which  is 
proposed  by  Zhang  and  Liu  in  [13]  [21].  The  electric  field 
can  be  expressed  as 

Ex  =  £%*  +  £%**,  (2) 

where  Einc  is  the  incident  field  from  an  infinite  line  current 
source  Jz 

EiJlc(p)  =  -jwn of  dqG[p,q)Jz{q).  (3) 

JD 


E*c at  is  the  scattered  field 

El^ipR.qr)  =  fD  dqG(pRy  q)A\{q)Ez  (q>qr)  (4) 

where  qr  and  pr  are  the  transmitter  and  receiver  positions. 
A k?{q)  =  k2(q)  -  and  k%  =  J*po€b  -  jupo^b-  D  is  the 
inhomogeneous  object  domain.  G  is  the  Greens  function 
for  a  homogeneous  background  medium 

G(p,q)  =  ±H?)(kb\p-q\).  (5) 

The  EBA  method  exploits  the  fact  that  G(p,q)  is  close  to 
5(p  —  q)  to  approximate  Eq.  4  as 

£|“*(Pfl>9r)  =  Ez(pR,qT)  [  dqG(j>R,q)Al(q).  (6) 

Jd 

Substituting  Eq.  6  into  Eq.  2,  we  obtain  the  solution  of  the 
electric  field 

Ex[p)  =  (1  -  f  dqG{p, q)Ae(q))-1Er-  (7) 

Jd 

The  approximation  of  the  scattered  field  is  obtained  by 
substituting  Eq.  7  into  Eq.  4.  With  the  improved  CGFFT 
method  ([13]  [21]),  the  approximate  scattered  electric  field 
can  be  computed  using 

Escat  = 

jW(p)  =  Ak?(l  -  T-l{T[0\T[Ak2]})  w 
B.  Reconstruction 

A  reconstructed  image  example,  shown  in  figure  lc,  is 
obtained  using  the  back-propagation  as  initial  solution 
followed  by  the  contrast  source  inversion(CSI  Berg[22]) 
method,  which  is  proposed  by  Zhang  and  Liu[13].  A  sum¬ 
mary  of  CSI  is  given  below.  First  define  the  contrast  source 
uj(j>)  as 

w(p)  =  x(p)Ez(p),  x(p)  =  Al(p)/k%.  (9) 

The  scattered  field  can  be  expressed  as 

£f“*(p)  =  kb  f  dpG(j>,q)x{q)Ez(q).  (10) 
Jd 

The  scattered  field  is  measured  in  the  domain  S  (outside 
D)  where  x  vanishes.  Assuming  no  measurement  error,  the 
measurement  data  r  satisfies 


The  contrast  source  can  be  obtained  by  minimizing  the  cost 
function  F(w,x)>  which  is  the  summation  of  the  error  in 
the  data  equation  and  state  equation. 

FM  =  (Ef=i  INI!)-1  Ef=i  \\rj  -  GsNI2+ 

(Ef=i  WxETjWlr1  Ef=i  UxEir  -  +  xGwiWl, 

(14) 

The  Polak  —  Ribiere  conjugate  gradient  procedure  is  used 
to  update  both  u  and  x  alternately. 

III.  Model  of  uncertain  permittivity  tissue 
IMAGE 

In  this  section  the  uncertainties  in  the  spatial  distribu¬ 
tion  of  the  breast  permittivity  properties  axe  modeled  sta¬ 
tistically.  First,  the  random  model  for  the  tissue  permittiv¬ 
ity  image  in  general  is  presented.  Then  the  random  model 
that  designates  both  the  normal  tissue  and  the  malignant 
tumor  is  presented. 

A.  Model  of  tissue  permittivity  image 

We  utilize  the  Markov  Random  Field  (MRF  [19])  to 
model  the  spatial  uncertainties  of  the  breast  permittivity 
image.  The  MRF  model  implements  the  idea  that  the  per¬ 
mittivity  of  spatially  nearby  tissue  is  similar.  It  does  not 
model  the  detailed  deterministic  anatomical  structure.  The 
Markov  property  captures  the  local  nature  of  the  permittiv¬ 
ity  distribution  and  thereby  reduces  the  uncertainty  space. 
A  Markov  Random  Field  satisfies  [23] 

p(xij\xki  €  ft  \  {(*,,?)})  =  p(xij\xki  € Nij)y  (15) 

where  Nij  is  the  neighborhood  system  of  (ij). 

There  are  several  MRF  models,  such  as  Auto-Models, 
Multi-Level  Logistic  Models.  The  Gaussian  MRF 
(GMRF,[24][23],  one  of  Auto-Models),  is  chosen  because 
it  can  easily  be  incorporated  into  the  Baysian  approaches. 
It  can  be  defined  as 


p(xij\xki  €  0\{(i,i)})  = 


where  pij  and  ofj  axe  the  mean  and  variance  at  (ij)-  and 
are  the  interaction  coefficients,  which  reflect  the 
local  characteristics.  The  matrix  form  of  P(ij)y(ki)  can  be 
expressed  as 


r(p)  =  $  f  dpG(pyq)x(q)Ez(q),P  €  S  (11) 
Jd 

We  have  the  state  equation  and  the  data  equation 

u(p)  =  xE?°  +  xGdv,p€  Dy  (12) 

r(p)  =  GsuiypeSy  v  ; 

where  G&u;  and  Gs&  are  defined  as 

GD,sv  =  k$fDdpG(p,q)x(q)Ez(q),  or  p£S 


A 0,0  A),  1  Po,M 

Pi,o  Ai,i  —  /Jim 

0M,O  0M,1  —  Am,M 

where  Aij,fci  has  been  simplified  to  assuming  a 

stationary  and  symmetric  property.  The  A  matrix  provides 
two  pieces  of  information:  one  is  the  neighborhood  area 
of  (ij),  which  is  a  (2 M  +  1)  x  (2 M  +  1)  square.  This 
implies  that  the  permittivity  at  (i  j)  only  depends  on  the 
permittivity  in  this  neighborhood  area.  Another  piece  of 
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information  is  the  influence  of  the  neighbors  on  the  point 
(ij),  as  represented  by  the  sign  and  value  of  the  interaction 
coefficients.  The  0  matrix  is  also  an  equivalence  to  the 
neighborhood  system  M  in  the  GMRF.  Using  conventional 
notation,  AT  in  this  model  is  composed  of  one-site  and  two- 
site  cliques  of  M-th  order.  A  clique  is  a  set  of  points  that 
are  neighbors  of  one  another. 

Eq.  16  only  defines  the  conditional  distribution  of  tissue 
permittivity.  According  to  the  Hammersley- Clifford  the¬ 
orem,  Eq.  16  is  equivalent  to  the  joint  distribution  of  a 
Gaussian  Field[25] ,  expressed  as 


p(X) 


X(k)  =  er(i,j),  k  =  (i  —  1)  x  N  +  j,l  <  i,j  <  N , 

(iS) 

where  X  is  a  N 2  x  1  vector,  an  re-ordered  mapping  of  per¬ 
mittivity  image  er  defined  on  an  N  x  N  lattice  fl.  The  sam¬ 
ple  space  of  €r  is  R1*2 .  \i  is  the  mean  vector  with  size  N2  x  1. 
Q  is  the  inverse  covariance  matrix  with  size  N2  x  N2.  Q 
should  be  symmetric  and  positive  definite. 

To  describe  the  equivalence,  the  relationship  of  the  pa¬ 
rameters  in  Eq.  16  and  those  in  Eq.  18  are  summarized. 


Eq.  18 
n(k) 

Q  =  UBtl 
n(fci,fc2) 


Eq .  16 

k=(i-l)xN  +  j 

kl  —  k2  =  (i  —  1)  *N  +  j 
*  Or  kl  £  k2 


B(kl,k2) 


s  0(ij),(ki)j  (M)  Z^fij 

{  0,  (kl)  gtfij 

kl  =  (i  —  l)xN+j 
k2=(k-l)x  N  + l. 


(19) 

When  the  random  field  is  assumed  stationary  and  isotropic, 
II  is  simplified  to  —J,  and  B  is  a  Toeplitz  block  circulant 
matrix. 


B.  Models  of  normal  and  malignant  tissue  permittivity  im¬ 
age 

Typically  in  many  studies,  malignant  and  normal  tissues 
have  been  modeled  as  homogeneous.  Here,  we  model  inho¬ 
mogeneous  tissue  characteristics  by  using  means,  variances 
and  interaction  coefficients  to  designate  both  the  normal 
tissue  and  the  malignant  tumor.  The  mean  values  of  tu¬ 
mor  and  background  are  fH  and  /if,  respectively,  which  axe 
scalar  variables.  The  interaction  coefficients  of  tumor  and 
background  are  matrix  0%  and  0b  respectively  with  assumed 
correlation  length.  An  algorithm  proposed  by  Rue  [26]  is 
used  to  fit  the  interaction  coefficients  to  the  Gaussian  field 
with  proposed  correlation  length. 

Transferring  the  local  characteristics  to  the  global  char¬ 
acteristics  for  our  Bayesian  signal  detection  approaches, 
we  have  a  uniform  mean  vector  pH°  =  {fibi  :  i  e  Q} 
for  Ho  condition,  and  a  non-uniform  mean  vector  jiHl  = 
(fH~fM>)xV  +  VHo  for  Hi  condition.  V  is  the  tumor  mean 


function  used  to  describe  the  position  S,  the  shape  and  the 
area  A  of  the  tumor.  We  consider  both  sharp  and  smooth 
tumor  mean  functions.  The  sharp  tumor  mean  function  is 
given  by 


V(r)  =  l,r  €  A(S) 
V(r)  =  0  ,r<£A(S) 


(20) 


The  smooth  tumor  mean  function  is  given  by 


V(r)  =  ^(r-S)4-^(r  — S)2  +  l,reA(S) 

D  =  n/0.6041, 

V(r)  =  0  ,r<£A(S) 

(21) 

where  S  is  the  tumor  position,  at  the  center  of  A.  n  is  a 
constant.  Both  the  sharp  and  smooth  function  are  assump¬ 
tions  since  real  data  of  tumor  spatial  distributions  is  not 
available  yet.  The  data  in  [5]  suggests  a  sharp  mean  func¬ 
tion,  i.e.  permittivity  at  the  tumor  boundary  that  is  higher 
than  that  at  the  tumor  center. 

From  the  local  characteristics  0b  and  variance  a 2,  we  de¬ 
rive  the  inverse  covariance  matrix  QH°  under  Hq  condition, 


=  %,),(«), 

=  0  ,(K)*AT« 


(22) 


Under  Hi  condition,  the  inverse  covariance  matrix  QHl, 
is  derived  from  both  0b  and  /3*,  as  well  as  cr2,  which  is  a 
function  of  tumor  position  S. 


QH'((ij),(kl))  - 

=  0,  (kl)  £  Mij 


(23) 

In  summary,  for  the  benign  tissue,  the  Ho  condition,  the 
distribution  of  er  can  be  expressed  as  er  ~  N(pcH°,QH°). 
For  the  malignant  tissue,  the  Hi  condition,  the  distribution 
of  er  can  be  expressed  as  er  ~  N  (fig1 ,  <3$ x).  I*1  the  follow¬ 
ing  sections,  we  suppress  the  superscript  HI  and  HO  if  it 
does  not  cause  any  ambiguity.  N  denotes  a  multivariate 
real  Gaussian  distribution. 

The  possible  parameter  space  is  limited  in  order  to  have  a 
definite  positive  covariance  matrix  of  the  Gaussian  Markov 
Random  field[27] .  However,  we  still  have  flexible  choices  for 
the  0%  <r’s  and  /i’s  in  order  to  capture  the  tissue  variance 
in  individuals  as  well  as  the  background  structure  variance 
across  patients.  The  model  is  simple  so  that  computational 
costs  are  tractable.  Since  the  Gaussian  field  is  a  class  of 


Gibbs  distribution  (Gibbs  [28]  [23]),  it  is  possible  to  extend 
the  GMRF  to  a  wide  range  of  Gibbs  distributions  to  model 
more  complex  microwave  tissue  characteristics. 


IV.  Optimal  detection  and  localization 

FRAMEWORK 

According  to  signal  detection  theory,  the  optimal  detec¬ 
tor  for  doubly  composite  hypotheses  is  the  likelihood  ratio 
(A)  of  the  data  vector  X  followed  by  a  threshold  whose  value 
is  determined  by  the  optimum  criterion  (T.  G.  Birdsail). 


A(Z)  = 


V(X  \  Hi) 
V(X  \  Ho) 


(24) 
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The  optimal  localization  processor  computes  the  a  posteri¬ 
ori  probability  of  the  tumor  location,  given  the  data  vector 
r  given  the  condition  frequently  we  use  the  maximum 
a  posteriori  probability  of  the  tumor  position  in  order  to 
obtain  quantitative  localization  performance  in  terms  of 
the  probability  of  correct  localization  (PCL). 

<25> 

Based  on  signal  detection  theory,  we  obtain  the  optimal 
physics  (tissue)-based  Bayesian  image  processor  using  two 
kinds  of  data:  the  data  of  the  defining  tissue  permittivity 
image  as  well  as  the  data  of  the  scattered  EM  field  mear 
surements.  The  Bayesian  results  are  presented  in  the  form 
of  the  likelihood  ratio  and  a  posteriori  probability  images. 
We  also  present  non-Bayesian  approaches  for  comparisons. 
For  all  these  approaches,  we  assume  in  this  paper  that  tu¬ 
mor  position  is  a  random  variable  with  uniform  distribution 
on  the  2D  lattice. 


Original  tissue  •  Iorv 

permittivity  image  ;  ~  - 


forward  propagation 


Or 


Reconstruction  of  tissue 
permittivity  image 


A)  Bayesian  Processor  for 
Uncertain  Permittivity  Image 

i)  Known  parameters 
(BPUPI) 

ii)  Estimated  parameters 
(EBPJJPI) 

B)  Threshold  Image  Processor 

(TIP)  _ 


'umor  detection  and  localization 


1  > 

C)  Bayesian  Processor  for 
Scattered  EM  field  of 
Uncertain  Pennitivity  Image 
(BPSEUPI) 

D)  Low  order  approximation  of 
BPSEUPI 
(LBP_SEUPI) 


umor  detection  and  localization 


Fig.  2.  Signal  detection  and  localization  using  uncertain  tissue  per¬ 
mittivity  models  and  EM  scattered  measurements 


Figure  2  shows  a  flow  graph  of  the  various  compar¬ 
isons.  The  left  part  in  figure  2  illustrates  the  first  Bayesian 
approach  for  tumor  detection  and  localization,  in  which 
the  reconstructed  image  is  optimally  post-processed.  The 
Bayesian  Processor  is  derived  using  the  mean,  variance  and 
interaction  coefficients  of  the  tissue  permittivity  image,  and 
its  output  is  called  the  Bayesian  Processor  for  Uncertain 
Permittivity  Image  (BP-UPI) .  It  is  assumed  that  the  re¬ 
construction  is  a  perfect  inverse  function  of  the  forward 
propagation,  i.e.  we  can  use  the  GMRF  model  presented 
in  sectionlll  for  both  original  and  reconstructed  permittiv¬ 
ity  image. 

The  performance  upper  bounds  are  obtained  based  on 
the  BP_UPI  algorithm  since  this  optimally  process  the  orig¬ 
inal  permittivity  image  data,  with  its  prescribed  uncertain¬ 
ties.  If  accurate  a  priori  knowledge  of  the  GMRF  is  not 
available,  a  suboptimal  approach  could  use  the  estimated 


parameters.  We  call  the  output  for  such  a  case  the  Esti¬ 
mated  Bayesian  Processor  for  Uncertain  Permittivity  Im¬ 
age,  EBPJJPI.  One  of  the  performance  lower  bounds  is 
obtained  by  defining  the  Threshold  Image  Processor(TTP) 
of  the  reconstructed  tissue  permittivity  image,  which  is  per¬ 
haps  similar  to  what  a  human  observer  might  do  in  looking 
at  such  an  image. 

A  second  approach  is  presented  in  the  right  part  of  Fig¬ 
ure  2.  Here,  the  likelihood  ratio  and  a  posteriori  prob¬ 
abilities  are  formed  directly  from  the  measurement  data, 
incorporating  the  a  priori  knowledge  of  the  uncertainties 
of  the  original  tissue  characteristics  as  projected  through 
the  nonlinear  propagation.  The  output  of  this  processor 
is  called  the  Bayesian  Processor  for  Scattered  EM  field  of 
Uncertain  Permittivity  Image(BP-SEMUPI).  The  Monte 
Carlo  Integration  is  used  to  make  the  computation  of 
BPJSEMUPI  tractable.  Since  the  scale  of  the  uncertainty  is 
very  large,  the  Monte  Carlo  Integration,  although  useful, 
takes  a  long  time  to  converge.  To  reduce  the  computa¬ 
tional  complexity,  we  propose  a  lower  order  approximation 
of  the  BP.SEMUPI,  called  the  LBPJ3EMUPI.  It  is  based 
on  projecting  the  uncertainties  from  the  measurement  do¬ 
main  onto  a  multivariate  Gaussian  field. 


A.  Bayesian  Processor  for  Uncertain  Permittivity  Image 
A.l  BP-UPI  detector 

In  the  problem  of  unknown  tumor  position  S,  the  likeli¬ 
hood  ratio  is 


A(«r) 


p(MHo) 


(26) 


The  a  priori  pdf  of  the  permittivity  image  has  been  mod- 


e!e<1  p(cr\HuS)  ~  N(jxs,Qs ) 

p(er\H0)  ~  N(ji,Q), 


(27) 


The  likelihood  ratio  for  the  BP-UPI 


A(«r)  «  Esen 


I  (*« - - - M.OTQ57(Ci - M«?) 

I  Qs  I  e  5 

(28) 


A. 2  BP-UPI  localization  processor 

The  a  posteriori  probability  of  the  tissue  permittivity, 
try  given  that  the  tumor  is  present  is  expressed  as 

Ps|£r(S|er)  S=  (29) 

Using  p(er  \  Hi,S)  ~  N(ps,Qs )  and  substituting  this  in¬ 
formation  into  Eq.  28  and  simplifying,  we  obtain 

p(S  |  er)  a  (er  -  p)T  Q{t r  -  m)+  /3Qn 

log(\  Qs  |)  -  («r  -  Ps)TQs(t r  ~  Vs)- 

With  some  careful  simplifications,  the  computations  will 
be  proportional  to  the  tumor  area  A. 

B.  Threshold  Image  Processor 

The  Threshold  Image  Processor  (TIP)  mimics  a  routine 
visual  examination  of  the  maximum  brightness  of  the  ob¬ 
ject  image. 
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B  J  TIP  detector 

The  threshold  detector  is  defined  by 

A  =  max  tr  (31) 

where  er  is  the  permittivity  image  data. 

B.2  TIP  localization  processor 

The  threshold  localization  processor  is  defined  by 

<32) 


C.2  BPJ3EMUPI  localization  processor 

In  the  measurement  domain,  we  derive  the  a  posteriori 
probability  of  the  tumor  position  S  given  the  received  mea¬ 
surements  r. 

P(5|r)  =  /£r  p(r/er,  S)p(er\S)der  (37) 

Substituting  Eq.  34  into  Eq.  37,  we  obtain 

p(S|r)  a  f€  e  2*«  p(€r|i7i,5)der. 

Monte  Carlo  integration  can  be  used  to  make  the  compu¬ 
tation  tractable 


where  S  is  the  tumor  position. 


C.  Bayesian  Processor  for  Scattered  EM  field  of  Uncertain 
Permittivity  Image 

C.l  BPJSEMUPI  detector 

In  the  measurement  domain,  The  likelihood  ratio  is  de¬ 
fined  by 


A(r) 


ES€n  /<„  P(r|Hi,er,S)P(£r|H1,S)p(S)dCr  /00, 
4  p(r\tr,HoMer\Ho)dtr  ’  W 


where  S  is  the  unknown  tumor  position. 

Incorporating  the  a  priori  knowledge  of  the  additive 
noise,  n  ~  calN(0,2o%I),  the  conditional  probability  den¬ 
sity  function  of  r  given  T(er  \  C)  is 


2  (r-T(erlC))H(r-T(eT.|C)) 

p(r  |  T(er),  C)  =  (2iralTK  e  (34) 


where  T  is  the  mapping  function  from  the  permittivity  im¬ 
age  to  the  measurement  data  and  T(er)  is  the  pure  mea¬ 
surement  signal  s. 

The  complete  likelihood  ratio  is  given  by 


p(Sk)oc]Te 


-33 


fc=  1 


(39) 


where  €rk(S)  follows  the  distribution  N(ps,  Qs ),  given  S.  L 
is  a  number  large  enough  to  make  the  Monte  Carlo  converge 
to  a  good  approximation  of  the  a  posteriori  probability. 

D.  Lower  order  approximation  of  BPJSEMUPI 

Considering  the  huge  time  cost  for  computing  BP-SEMUPI, 
a  lower  order  distribution  is  considered  to  approximate  the 
true  distribution  of  the  measurement  data,  and  the  resul¬ 
tant  processor  is  called  LBPJ3EMUPI.  The  multivariate 
Gaussian  model  is  chosen  which  demonstrates  good  per¬ 
formance  in  the  simulations.  From  a  generated  large  data¬ 
base,  we  estimate  the  mean  and  covariance  matrix  of  the 
measurement  data  under  both  Hq  and  Hi,  which  are  M 
and  D  for  ff0  and  M$  and  D$  for  (i?i,S).  The  size  of  the 
Ms  and  Ds  matrices  are  the  same  as  the  size  of  the  possi¬ 
ble  tumor  positions.  These  estimates  are  then  used  in  the 
LBPJ3EMUPI  processors  to  make  detection  and  localiza¬ 
tion  decisions. 


D.l  LBP.SEMUPI  detector 

Using  a  multivariate  complex  Gaussian  model,  we  have 


A(r)°c 


Sse  n  fer 


o— TfcOV^r-ru,.)) 

*  a 


p(er  |  HuS)de r 


L 


(r-T(<r))H  (r-TM) 


p(€r  |  H0)de r 


(35) 

Because  T(er)  is  a  complex  nonlinear  function,  we  can  not 
obtain  a  closed  form  for  the  likelihood  ratio.  However, 
Monte  Carlo  Integration  can  be  used  to  compute  A(r)  in 
a  computationally  efficient  manner.  This  can  be  expressed 
as 


A(r)  oc 

EM  * 

3=1  ,Sjen  e _ _ 

EzZS  1  ))  " 

k  1  Hq  )  e  n 

(36) 

where  M  is  the  number  of  all  possible  tumor  positions  and 
L  is  a  number  large  enough  to  make  the  Monte  Carlo  Inte¬ 
gration  a  good  approximation  to  the  likelihood  ratio. 


p(Ttr\C))  = 


I  £c  1  ~(T.. -Mc)hDc(T.t  -Me) 

(27r)^C 


(40) 


where  subscript  C  represents  (Hq)  or  (Hi  ,S).  K  is  the  num¬ 
ber  of  sensors. 

Using  the  theorem  of  Multiplication, 


p(r  \C)=  [  p(r\  TCrlc)p(TcAc)dTtr.  (41) 


Substituting  Eq.  40  and  Eq.  34  into  Eq.  41,  we  obtain 
p(r  |  C)  =  ^ D° 1 ^Dci jr 1  cMgl  Dci  Mci  ~ ~M* Dc Mc , 


Dei  =2^  +  Z>c, 

Mci  =Dc{(£t  +  DcMc). 

(42) 
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The  likelihood  ratio  can  be  expressed  as 


wx  ZszoPir  I  H1,S)p(S) 

X(T)  -  wm 


(43) 


where  M  is  the  number  of  all  possible  tumor  positions. 

Substituting  Eq.  42  into  Eq.  43,  with  C  replaced  by 
(#i,S)  and  {H0)  respectively,  we  obtain  the  simplified  form 
of  the  LBP  J3EMUPI  detector 


A(r)  °c  Esc n  I  Ds  II  Dsi  I  1 

gMsi**  DsiMsi — MsH  Ds  Ms  —MiH  D1M1 


(44) 


where,  Msu  Dsi,  M\  and  D\  follow  Eq.  42. 

When  approaches  zero,  i.e.  there  is  no  measurement 
error,  the  expressions  are  simplified  to 

p(r  |  C)  =  J^J.e-(r-Mc)HDc(r-Mc), 

A(r)  OC  Y,S€Cl  I  As  |*  e(r-M)*D(r-M)-(T-Ms)BDs(T-Ms) 

(45) 

D.2  LBPJ3EMUPI  localization  processor 


and  background  are  /3 *  and  to  fit  the  Gaussian  field  with 

3  and  30  pixel  correlation  lengths. 


1.0000 

-0.1879 

0.0191  ‘ 

1  0.1  r 

0.1  0.1 

-0.1879 

-0.1724 

0.0453 

A  = 

0.0191 

0.0453 

0.0007 

(48) 


We  also  consider  the  situation  where  both  tumor  and  nor¬ 
mal  tissue  have  the  same  interaction  coefficient  /3b  in  the 
first  simulation  scenario.  Two  tumor  types,  one  with  sharp 
boundaries  and  the  other  with  smooth  boundaries  mean 
are  used  in  the  first  scenario,  as  shown  in  figure  3a  4a. 

B.  Performance  evaluation 

To  evaluate  the  detection  performance,  we  use  the  Re¬ 
ceiver  Operating  Characteristic  (ROC),  which  is  a  plot  of 
the  Probability  of  Detection  {Pd)  as  a  function  of  Proba¬ 
bility  of  False  alarm  (Pf)*  Pd  and  Pp  are  obtained  from 
the  probability  density  functions  of  the  likelihood  ratio  A 
under  each  hypothesis  in  that 

Pd  =  XT p(A  I  Hi)d\  PF  =  p(A  I  Ho)dX  (49) 


The  localization  problem  has  already  been  defined  in  Eq. 
37.  With  the  multivariate  Gaussian  model  approximation, 
the  conditional  probability  density  function  of  r  given  Hi 
and  S  is  defined  in  Eq.  42,  with  C  replaced  by  (ifi,S).  Then 
the  a  posteriori  probability  expression  is  given  by 

p(S  I  r)  «  Nye<Ps,%'^~Mj,PsMf, 

Dsi  =  ^r  +  Ds,  (46) 

Ms i  =  Dsi{^r  +  DsMs)- 

When  there  is  no  additive  noise,  the  LBPJSEMUPI  local¬ 
ization  processor  is  simplified  to 

p(£  |  r)  oc  log{\  Ds  |)  —  (t-  —  Ms)ffDs(r  -  Ms).  (47) 

V.  Simulation  Results 
A .  Basic  parameters 

In  the  simulations  in  this  paper,  the  tissue  permittivity 
image  has  size  9.2cm  x  9.2cm,  529  pixels  and  the  unit  area 
is  4mm  x 4mm.  For  the  first  three  simulation  scenarios, 
there  are  24  sensors  evenly  spaced  outside  the  object  area, 
7  at  each  side  and  1.6cm  between  the  neighborhood  sensors. 
The  frequency  used  for  propagation  and  reconstruction  is 
1GHz.  The  last  simulation  scenario  discusses  the  effect  of 
different  array  configurations. 

The  mean  value  of  normal  background  tissue  was  set 
to  fib  —  30,  which  is  within  the  range  of  values  reported 
experimentally  for  real  breast  permittivity  [6].  The  mean 
value  of  malignant  tumor  /j *  is  considered  to  be  a  variable 
in  the  first  simulation  in  order  to  study  the  effect  of  the 
tumor  contrast,  and  it  takes  constant  value  of  40  in  the 
later  simulations.  The  variance  for  both  the  background 
and  the  tumor  is  1.  The  interaction  coefficients  of  tumor 


We  also  use  the  detectability  index,  d,  sometimes  called 
the  deflection  coefficient,  to  evaluate  the  detection  perfor¬ 
mance.  The  scalar  variable  d  completely  characterizes  the 
ROC  for  the  case  where  the  underlying  probability  density 
functions  are  Gaussian.  In  other  cases,  d  can  be  estimated 
from  the  ROC  curve  through  Eq.  50  using  a  data  fitting 
method  and  used  as  an  approximation  to  the  ROC’s. 

PD  =  1  -  Gau2cdf{Gau2inv{l  —  Pp)  —  d)  (50) 

where  Gau2cdf  and  Gau2inv  are  normalized  Gaussian  cu¬ 
mulative  density  function  (CDF)  and  inverse  CDF  func¬ 
tions. 

To  evaluate  the  localization  performance,  we  use  Proba¬ 
bility  of  Correct  Localization  (PCL).  PCL  is  the  ratio  of  the 
number  of  correct  localization  divided  by  the  total  number 
of  trials. 

A  database  of  the  original  and  reconstructed  tissue  per¬ 
mittivity  image  and  measurement  data  is  generated  for  the 
performance  evaluation,  including  400  realizations  with  a 
tumor  at  random  positions  and  400  realizations  with  no  tu¬ 
mor.  The  test  database  is  independent  of  a  large  database 
generated  to  study  the  statistics  of  the  measurement  data. 

C.  Performance  upper  bound  and  tumor  characteristics 

The  detection  performances  of  BP_UPI  are  shown  in  fig¬ 
ures  3b  3c,  4b,  and  4c.  In  those  figures  the  detectabil¬ 
ity  is  plotted  as  a  function  of  tumor  area,  tumor  con¬ 
trast  ^hi-j+ho  ?  sharp  or  smooth  tumor  mean  function,  and 
whether  or  not  the  interaction  coefficients  for  malignant 
and  background  tissue  are  the  same  or  different.  The  com¬ 
putation  follows  Eq.  28. 

The  localization  performances  are  shown  by  PCL  curves 
in  figures  3d,  3e,  4d,  and  4e.  That  computation  follows  Eq. 
30. 

The  simulation  results  demonstrate  that 
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1. ^The  tumor  detection  and  localization  performance  in¬ 

s'  1  HO 

creases  with  the  tumor  contrast  *= — ^ — ,  which  is  shown 
in  all  8  plots. 

2.  By  comparing  figure  3  and  figure  4,  it  can  be  seen  that 
it  is  easier  to  detect  and  locate  the  tumor  with  a  sharp 
weight  function  than  that  with  a  smooth  weight  function. 

3.  By  comparing  (b)  and  (c),  (d)  and  (e)  in  both  figures 
3  and  4,  one  can  see  that  capturing  different  interaction 
coefficients  for  the  background  and  the  tumor  will  help  to 
improve  the  detection  and  localization  performance. 

4.  The  performance  as  a  function  of  tumor  area  is  not  a 
simple  scaling  relationship  as  that  of  tumor  contrast  since 
it  is  influenced  by  the  tumor  mean  function  and  the  in¬ 
teraction  coefficients  together.  When  the  tumor  and  the 
background  have  the  same  interaction  coefficients  and  the 
tumor  has  a  smooth  shape  function,  the  larger  tumor  area 
may  cause  stronger  similarity  between  the  tumor  and  the 
background  such  that  the  detection  performance  may  de¬ 
crease,  as  shown  in  4c, e). 


Fig.  3.  Detection  and  localization  performance  of  BP.UPI  using 
original  tissue  permittivity  image  as  a  function  of  tumor  contrast, 
size,  local  characteristics,  with  sharp  tumor  mean  function.  a)an 
example  of  sharp  tumor  mean  function,  tumor  at  (11,7),  occupy¬ 
ing  25  pixels.  b)c)ROC,  d)e)PCL  curves.  Interaction  coefficients 
for  background  and  tumor:  left-differ, right-same.  From  top  to  bot¬ 
tom, contrast  Left- [5  2.5  2  1.5  1  0],  Rigbt-[5  2.5  2  1.5  1] 

D.  Convergence  of  the  covariance  matrix  estimation 

The  estimation  of  the  mean  and  covariance  matrix  is  an 
important  part  for  the  performance  of  the  computation¬ 
ally  simpler  EBPJJPI  and  LBPJSEMUPI.  Since  the  un¬ 
certainty  scale  is  large,  we  need  a  large  number  of  samples 
for  the  estimates  to  converge. 


Fig.  4.  Detection  and  localization  performance  of  BP-UPI  using 
original  tissue  permittivity  image  as  a  function  of  tumor  contrast, 
size,  local  characteristics,  with  smooth  tumor  mean  function.  a)an 
example  of  smooth  tumor  mean  function,  tumor  at  (11,17), occupy¬ 
ing  25  pixels.  b)c)ROC,  d)e)PCL  curves.  Interaction  coefficients 
for  background  and  tumor:  left-differ , right-same .  From  top  to  bot¬ 
tom, contrast  PKl-PgQ  :  Left-[5  2.5  2  1.5  1  0],  Right-[5  2.5  2  1.5  1] 

Two  approaches  are  used  to  study  the  convergence  of  the 
estimates, 

1.  Determine  the  performance  of  the  processor  versus  the 
number  of  samples  used  to  form  the  estimates.  The  ROC 
curve  should  reach  some  stable  state  as  the  number  of  the 
samples  increases. 

2.  Estimate  the  mean  and  covariance  matrix  of  the  origi¬ 
nal  tissue  permittivity  image  data  from  the  database,  and 
compare  them  with  the  true  values.  The  difference  should 
be  within  a  small  threshold.  The  number  of  unknowns  is 
529  for  the  tissue  permittivity  image  and  600  for  the  mea¬ 
surement  data.  Since  they  are  comparable,  we  can  use  the 
convergence  of  estimates  using  the  image  data  to  predict 
that  using  the  measurement  data. 

Figure  5  illustrates  the  detection  performance  of  the 
BP_UPI  and  LBP_SEMUPI  detectors  versus  the  number 
of  samples  used  in  the  estimations  of  mean  vector  and  co- 
variance  matrix.  Figure  5a  uses  the  original  permittivity 
image  data  and  it  can  be  seen  that  the  ROC  curves  be¬ 
come  stable  near  the  perfect  upper  left  comer.  The  corre¬ 
sponding  deflection  coefficient  converges  to  around  6.4  af¬ 
ter  1400  samples.  It  should  be  noted  that  6.4  is  the  upper 
limit  of  the  deflection  coefficient  that  we  cam  estimate  using 
the  data  fitting  method.  Figure  5b  uses  the  scattered  EM 
measurement  data.  The  ROC  curves  converge  with  some 
degradation  compared  with  those  using  the  original  permit¬ 
tivity  image  data.  The  corresponding  deflection  coefficient 
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converges  to  around  3.6  after  1400  samples.  The  perfor¬ 
mance  degradation  as  sample  size  decreases  also  reflects 
the  effect  of  non-accurate  a  priori  knowledge  of  the  mean 
and  covariance  matrix  on  the  processor’s  performance. 


Fig.  5.  The  convergence  of  the  detection  performance  as  a  function 
of  the  number  of  samples  used  to  estimate  the  mean  and  covariance 
matrix.  Assume  the  tumor  position  unknown,  a)  ROC  curves  of 
BP-UPI  detector  using  the  original  tissue  permittivity  image  data  b) 
ROC  curves  of  LBP-SEMUPI  detector  using  the  measurement  data 

Table  I  lists  the  value  of  the  estimated  coefficients  from 
2000  samples  and  the  true  values  using  the  original  tissue 
permittivity  image  data.  There  are  35%-50%  increases  in 
the  absolute  estimation  values.  However,  the  performance 
curves  demonstrate  that  this  level  of  accuracy  in  interaction 
coefficients  is  OK  for  a  stable  performance. 

TABLE  I 

The  true  and  the  estimated  interaction  coefficients 


A>n 

fibl2 

A>13 

fib  22 

ture 

1 

-0.1879 

0.0191 

-0.1724 

1.3557 

-0.2547 

0.0255 

-0.2337 

A>23 

Pb33 

wsm 

$12.22 

ture 

0.0453 

0.00067 

1 

0.1 

estimate 

0.0607 

0.0007 

1.3613 

0.1459 

E.  Performance  of  proposed  Bayesian  approaches 
E.l  An  example 

Figure  6a  is  an  example  of  a  stochastic  background  per¬ 
mittivity  image  of  the  tissue,  along  with  a  simulated  tumor, 
modeled  by  the  GMRF.  Figures  6b,  c,  and  d  are  the  recon¬ 
structed  permittivity  images  from  the  perfect  measurement 
data  or  from  60dB  and  50dB  noisy  measurement  data. 

The  signal  detection  approach  using  the  permittivity 
data  is  to  compute  the  a  posteriori  probability  of  the  tu¬ 
mor  location  given  either  the  original  permittivity  image, 
or  the  reconstructed  image  as  data,  with  Eq.  30.  Figure  6e 
gives  an  upper  bound  on  tumor  localization  by  plotting  the 
a  posteriori  probability  of  tumor  location  using  the  tissue 
data  of  Figure  6a.  Figure  6i  gives  a  sub-optimal  perfor¬ 
mance  comparison  for  tumor  localization  by  plotting  the 
a  posteriori  probability  of  tumor  location  using  the  tissue 
data  of  Figure  6a  and  the  estimation  of  the  mean  and  co- 
variance  matrix.  Figures  6f-h  shows  the  a  posteriori  plot  of 
tumor  location  based  on  post  processing  the  reconstructed 
tissue  data  shown  in  Figures  6b-d. 


Using  the  measurement  data,  the  a  posteriori  probability 
of  tumor  location  is  computed  with  Eq.  46.  Figure  6j-l 
shows  the  a  posteriori  plot  of  tumor  location  base  on  the 
same  measurements  to  get  the  reconstruction  in  figures  6b- 
d. 

Plots  6g-h  show  that  at  60dB  and  50dB  additive  noise 
condition,  the  BPJJPI  using  the  reconstruction  tissue  data 
misses  the  correct  tumor  localization,  and  figures  6k-l  show 
that  the  BPJ3EMUPI  using  the  measurement  data  gets  the 
correct  location  of  the  tumor  with  high  probability.  This 
is  a  specific  example  where  the  BP_SEMUPI  works  better. 
In  the  following  sections,  it  is  demonstrated  statistically 
that  the  processors  based  directly  on  the  likelihood  ratio 
and  a  posteriori  probabilities  of  the  measurement  data  have 
better  performance. 


Fig.  6.  a)  Original  tissue  permittivity  image;  b)-d)  Reconstructed 
permittivity  from  the  measurement  data  b)without  additive  noise 
c)with  additive  noise,  SNR=60dB  d)with  additive  noise,  SNR=50dB; 
e)-h)  the  a  posteriori  probability  of  the  tumor  position  given  the  per* 
mittivity  image  data  and  HI  condition-p(S  |  €r,  HI),  data  er  comes 
from  a)-d);  i)the  p(S  |  €r,  HI)  computed  using  the  estimation  of  the 
mean  and  covariance  matrix  and  the  data  from  a);  j)-l)  the  o  poste¬ 
riori  probability  of  the  tumor  position  given  the  measurement  data 
and  HI  condition-p(S  |  r,  Hi) 


E.2  Detection  performance 

Figure  7  illustrates  detection  performance  comparisons 
assuming  no  additive  noise.  In  figure  7  the  ROC  of  the 
threshold  detector  provides  a  performance  lower  bound. 
TIP  mimics  the  way  a  routine  visual  examination  of  the 
image  might  be  done.  Although  it  is  not  sophisticated, 
it’s  ROC  reflects  the  problem  of  high  positive  predictive 
Value(PPV)  of  conventional  mammography.  Using  the  re¬ 
constructed  permittivity  image  data,  the  BP-UPI  detector 
is  much  better  than  the  TIP  detector,  especially  when  the 
probability  of  false  alarm  is  low,  because  the  BP-UPI  uti- 
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lizes  the  a  priori  knowledge  of  the  tissue  background  across 
0  patients  and  the  a  prior  knowledge  of  the  different  charac¬ 
teristics  of  the  normal  and  malignant  tissues  to  improve 
the  detection  performance.  With  a  shared  large  database, 
accurate  a  prior  knowledge  will  be  available.  Even  with 
uncertain  a  priori  knowledge,  as  shown  in  figure  5,  the 
performance  is  still  satisfactory  compared  with  that  of  the 
threshold  method. 

Figure  7  also  shows  that  the  BP-SEMUPI  detector  us¬ 
ing  the  measurement  data  is  better  than  BPJJPI  detector 
using  the  reconstructed  permittivity  data,  yet  worse  than 
BPJJPI  detector  using  the  original  permittivity  data.  It 
should  be  noted  that  in  reality,  we  can  not  get  the  original 
tissue  permittivity  image  directly  but  the  scattered  EM 
measurements.  However,  it  provides  an  upper  bound  for 
performance  evaluations.  It  is  better  than  the  BP-SEMUPI 
because  the  forward  EM  field  maps  the  variables  from  the 
original  permittivity  domain  to  the  measurement  domain, 
which  reduces  the  random  variable  space  and  decreases  the 
detectability.  Proofs  in  the  Appendix  support  this  explana¬ 
tion.  In  addition,  it  is  because  of  the  multivariate  Gaussian 
approximation  used  to  compute  Eq.  37.  The  LBP_SEMUPI 
is  better  than  the  BPJJPI  using  reconstructed  permittivity 
data  because  of  the  inherent  algorithm  limitations  in  the 
reconstruction  procedure. 


Fig.  7.  Detection  performance  of  BPJJPI  using  original  tissue  data, 
BPJJPI  using  reconstructed  permittivity  image  data,  TIP  using  re¬ 
constructed  permittivity  image  data  and  BP-SEMUPI  using  mea¬ 
surement  data.  No  additive  noise  presents  at  sensors,  tumor  position 
unknown 


Fig.  8.  Detection  performance  comparisons.  Tumor  position  un¬ 
known.  SNft=(Inf,  60,50]  dB  a)  BPJJPI  and  TIP  detectors  using  the 
reconstructed  tissue  permittivity  image  b)  BPJJPI  detector  using  the 
reconstructed  tissue  permittivity  image  and  LBP-SEMUPI  detector 
using  the  measurement  data 

Figure  8  displays  the  performance  of  the  detectors  when 


noise  is  present.  Additive  noise  degrades  the  performance 
of  all  three  detectors.  Figure  8a  shows  that  for  the  three 
SNR  conditions,  the  BPJJPI  detector  is  better  than  the 
TIP  detector.  The  ROC  curve  at  60dB  SNR  condition 
is  so  close  to  the  ROC  curve  for  a  very  large  SNR  that 
we  consider  60dB  as  a  threshold  for  the  BPJJPI  detector 
using  reconstruction  data.  Figure  8b)  shows  that  for  the 
three  SNR  conditions,  the  LBPJSEMUPI  detector  using 
the  measurement  data  is  better  than  the  BPJJPI  detector 
using  the  reconstructed  tissue  permittivity  image  data. 

E.3  Localization  performance 

Figure  9  displays  the  localization  performance  of  the 
LBP-SEMUPI  using  the  measurement  data,  and  the 
BPJJPI  and  the  TIP  using  the  reconstructed  tissue  per¬ 
mittivity  data.  The  localization  performance  is  shown  by 

a)  window  diameter  =0  b)  window  cSameter  =1 


Fig.  9.  Localization  performance:  solid  line,  LBP-SEMUPI  using 
the  measurement  data;  dotted  line,  BPJJPI  using  the  reconstructed 
tissue  permittivity  data;  dashed  line,  TIP  using  the  reconstructed 
tissue  permittivity  data 

the  probability  of  correct  localization  (PCL)  curves.  The 
localization  is  correct  if  the  picked  position  is  located  within 
the  test  window.  If  it  is  required  that  the  picked  position 
is  exactly  the  same  as  the  real  position  to  be  true,  the  win¬ 
dow  diameter  is  zero.  For  the  other  values  of  the  window 
diameter:!, 2, 2.8284,  the  window  sizes  axe  9,13,25.  Figure  9 
demonstrates  that  at  50dB  or  60dB  SNR,  both  the  BPJJPI 
and  the  LBP_SEMUPI  have  performance  better  than  0.9 
PCL.  They  are  much  better  than  the  threshold  method 
for  all  the  tested  SNR  conditions.  The  LBP-SEMUPI  lo¬ 
calization  using  the  measurement  data  is  the  best  of  the 
three.  The  results  of  the  detection  and  localization  per¬ 
formances  show  that  the  LBP-SEMUPI,  which  processes 
the  measurement  data  directly,  provides  a  potential  way  to 
help  a  doctor  make  a  better  diagnosis. 

F.  Performance  and  Array  configurations 

This  simulation  scenario  studies  the  detection  perfor¬ 
mance  as  a  function  of  the  sensor  array  configuration.  Fig- 
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use  10  displays  the  deflection  coefficient  indexed  by  2D 
source  positions,  given  the  possible  tumor  position  is  known 
to  the  detectors. 


Fig.  10.  Deflection  coefficients  as  a  function  of  the  tumor  positions, 
assume  tumor  position  known  for  the  detector,  a)  BPJUPI  using  orig¬ 
inal  permittivity  data  b)  EBP-UPI  using  original  permittivity  data 
c)-f)  LBP.SEMUPI  using  measurement  data.  ^Configuration  1:1:24, 
24  sensors  work  d)  Configuration  1:2:24,  12  sensors  work  e)  Configu¬ 
ration  2:2:24,  12  sensor  work  f)  Configuration  1:1:12, 12  sensors  work 


Figure  10a  provides  the  upper  bound  by  computing  the 
deflection  coefficient  of  the  BPJUPI  using  the  original  per¬ 
mittivity  data.  Figure  10b  is  a  sub-optimal  comparison 
by  computing  the  deflection  coefficient  of  EBP  JCJPI  using 
the  original  permittivity  data.  The  degradation  in  figure 
10b  compared  with  that  in  figure  10a  reflects  the  effect  of 
the  estimation  procedure,  which  is  independent  of  the  tu¬ 
mor  positions.  Figures  10c-  lOf  is  the  performance  of  the 
LBP_SEMUPI  detector  for  four  different  array  configura¬ 
tions.  The  four  plots  show  a  strong  relationship  between 
the  detection  performance  and  the  distance  of  the  tumor 
to  the  sensors.  The  deeper  the  tumor  is  located,  the  more 
difficult  it  is  to  detect. 


Fig.  11.  Detection  performance  using  the  measurement  data  and 
the  reconstructed  permittivity  image  data.  Assume  tumor  position 
unknown.  Comparison  of  different  array  configurations. 

Figure  11  displays  the  ROC  detection  performance  when 
the  possible  tumor  position  is  unknown.  By  comparing  the 
ROC  curves  of  the  first  configurations  with  the  three  other 
configurations  in  figure  11,  one  sees  that  the  performance 
degrades  due  to  the  limited  number  of  sensors.  The  ratio 
of  the  data  length  to  the  uncertainty  variable  number  is 
600:529  for  the  first  configuration  and  156:529  for  three 
other  configurations.  On  the  other  hand,  if  we  lower  the 


requirement  on  the  processor’s  performance,  we  can  reduce 
the  necessary  number  of  sensors.  The  plots  also  show  how 
the  configuration  of  sensors  will  influence  the  processor’s 
performance.  For  the  case  where  the  tumor  position  follows 
a  uniform  a  priori  distribution,  it  appears  that  the  more 
evenly  spaced  the  sensors,  the  better  the  performance.  But 
if  the  a  priori  tumor  position  is  not  a  uniform  distribution, 
we  could  improve  the  performance  by  having  more  sensors 
close  to  the  possible  tumor  positions. 

VI.  CONCLUSION 

This  paper  considers  two  Bayesian  approaches  for  breast 
tumor  diagnosis,  one  post  processes  the  reconstructed  tis¬ 
sue  permittivity  image  data.  The  other  incorporates  the 
scattered  EM  measurement  data  into  an  optimal  likelihood 
ratio  detector.  The  breast  permittivity  cross  section  propa¬ 
gation  medium  is  modeled  by  a  Gaussian  Markov  Random 
Field.  The  simulations  compare  these  two  approaches.  The 
simulation  results  also  provide  an  upper  bound  on  early  de¬ 
tection  performance  and  localization  of  malignant  tissue, 
as  a  function  of  uncertain  and  variable  tissue  permittiv¬ 
ity  characteristics,  tumor  contrast,  tumor  size,  and  local 
characteristics  and  shape.  The  effectiveness  of  using  the 
multivariate  Gaussian  distribution  approximation  for  the 
measurement  data  to  reduce  the  computational  complexity 
was  also  illustrated.  Finally,  the  effect  of  array  configura¬ 
tions  on  detection  performance  was  illustrated. 
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Appendix:  Deflection  Coefficient  Inequalities 

For  linear  projection  and  reconstruction  medical  system, 
the  performance  of  the  optimal  detectors  to  distinguish  the 
tumor  of  some  fixed  position  from  the  normal  tissue  back¬ 
ground  can  be  precisely  represented  by  the  deflection  coeffi¬ 
cient  if  multivariate  normal  distributions  are  used  to  model 
the  object  image  and  the  measurement  error.  In  this  part, 
we  derive  the  deflection  coefficient  based  on  these  assump¬ 
tions  using  the  original,  measurement  and  reconstruction 
data,  and  prove  relationships  with  a  matrix  inequality  the¬ 
orem.  Furthermore,  the  effects  of  the  forward  propagation, 
the  reconstruction  procedure  and  the  measurement  error 
on  the  detection  performance  are  clarified. 

Theorem  1:  uTC~lu  >  uTBT(BCBT)~1Bu,  if  C  is  a 


NxN  real  symmetric  positive  definite  matrix,  u  is  a  Nxl 
vector,  B  is  a  MxN  matrix  and  BCBT  is  non  singular. 

Proof:  We  have  the  matrix  version  of  the  well-known 
Cauchy-Schwarz  inequality[29] 

XTX>LXTY{YTY)~lYTX ,  (51) 

where,  X  is  a  nxm  matrix,  Y  is  a  nxq  matrix,  and  YTY 
is  nonsingular.  >i  is  the  Loewner  ordering  relationship. 

As  C  is  a  real  symmetric  positive  definite  matrix,  we 
have  C  =  LLT  and  C~l  =  L~TL~l ,  where  L  is  a  NxN 
simple  matrix.  Substitute  X  —  L~T  and  Y  =  LTBT  into 
Eq.  51,  we  get  i.e.  uTC^u  > 

utBt(BCBt)'1Bu.  M 

Deflection  Coefficient  after  linear  transformation 

Any  linear  transformation  can  be  represented  by  a  ma¬ 
trix  A.  Assume  x  is  a  Nx  1  vector  of  multivariate  Gaussian 
distribution,  y  is  a  Mx  1  vector  transformed  from  x  by  MxN 
matrix  A 

y- Ax.  (52) 

Assume  the  binary  hypotheses  for  x  are 

HI:  x~N(0uO 

HO:  x~N{60,C),  [06) 

where  Oi  is  the  Nxl  mean  vector;  C  is  the  NxN  covariance 
matrix,  which  is  real,  symmetric  and  positive  definite. 
From  Eq.  52  and  Eq.  53,  we  have 

E{y  |  FI}  =  Ah 

E{y\HQ}  =  A00  (54) 

COV{y}  =  ACAt. 

We  derive  the  deflection  coefficient  no  matter  whether  the 
COV(y)  be  singular  or  non-singular.  Let  Ra  =  rank  (A). 
There  are  Ra  independent  rows  in  A.  We  assume  Ai  is 
composed  of  Ra  independent  rows  of  A,  and  A2  is  com¬ 
posed  of  the  rest  of  the  rows  of  A.  A 2  can  be  represent  by 
a  linear  combination  of  rows  in  Ax,  i.e.  A?  =  GA\.  Then, 
we  have 


2/2  =  A2X  =  GA\x  —  Gy\. 


From  Eq.  55,  we  get  py(y)  =  S(y2  -  Gyi)pVl(yi).  Since  y2 
is  determined  by  t/i,  which  is  independent  of  the  hypothe¬ 
ses,  we  find  that  the  likelihood  ratio  of  y  is  the  same  as 
the  likelihood  of  yi ,  and  so  is  the  deflection  coefficient.  It 
can  be  proved  that  GOV  (yi)  =  AiCAj  is  a  real  symmet¬ 
ric  positive  definite  matrix.  Then  the  post-transformation 
deflection  coefficient  becomes 

d?  =  A0TA?  {A1CA^)~1A1A0.  (56) 

Deflection  coefficients  using  three  kinds  of  data 
1.  Original  image  data 

The  binary  hypotheses  have  been  defined  in  Eq.  53,  where 
HI  represents  the  normal  tissue  background  with  a  tumor 
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present,  HO  represents  the  normal  tissue  background  only, 
and  x  in  this  problem  is  the  permittivity  data  er  used  for 
detection.  The  deflection  coefficient  with  the  original  ob¬ 
ject  data  becomes 


A  0T  =  Oi  —  6q. 


(57) 


2.  Measurement  data 

If  we  consider  the  measurement  error,  and  use  the  additive 
independent  Gaussian  vector  to  represent  it,  the  measure¬ 
ment  data  r  can  be  represented  by 

r  =  Aer  +  n  =  [  A  IM  ]  £  J  ]  (53) 

where  A  is  a  MxN  matrix,  /^isaMxM  Identity  matrix, 
er  is  an  Nxl  vector,  n  is  a  Mxl  vector.  Using  Eq.  55,56, 

with  A  replaced  by  [A  Im],  x  replaced  by  |  j ,  and  y 

replaced  by  r,  we  obtain  the  deflection  coefficient  using  the 
measurement  data 


d?m  =  MTA[(A1CAj  AO,  (59) 


where,  [Ai  Imx]  are  composed  of  Ra  independent  rows  of 
matrix  [A  Im],  Ra  =  rank([A  Im])- 
3.  Reconstructed  image  data 

After  reconstruction,  the  estimated  object  data  €T  can  be 
represented  by 


Cr  =  Br  =  B 


n 

1*2 


-  Dru  D  = 


B 

BG 


(60) 


where,  n,r2  correspond  to  yi,j/2  in  Eq.  55.  G  is  a  ma¬ 
trix  that  satisfies  r2  =  Gr\,  Using  Eq.  [55,56]  again,  with 
A  replaced  by  D,  x  replaced  by  r\  and  y  replaced  by 
we  obtain  the  deflection  coefficient  of  the  optimal  detector 
using  reconstruction  data 


4 


A0t(B1A1)t(B1A1C(B1A1)t 
+B1Bfo*)-1B1A1  A6, 


where  B\  is  composed  of  Rb  independent  rows  of  D,  Rb  = 
rank(D). 


Inequality  relationships  of  deflection  coefficients 

We  have  d%  —  A0TC~1A6.  Assume  0%  —  0,  we  have 
=  &0TA^ (AiCAf  )“*  Ai  AO.  Apply  theorem  1,  we  get 
dP0  >  It  can  be  proved  that  is  inversely  propor¬ 
tional  to  cr%,  or  d^  <  d^,  so  we  get 

«€>£.  (62) 

Let  A©  =  A\A6  and  E  =  AiCA{ and  substitute 
them  into  Eq.  59  and  Eq.  61,  we  have 


a  eTE~1Ae 

AeT(B1)T{B1EBf)~1BiAe. 


(63) 


Since  E  is  real  symmetric  positive  definite,  BEBT  is  non 
singular,  apply  theorem  1,  we  get 

d2m>4.  (64) 


Two  special  cases 

1.  Case  1:  Assume  the  measurement  error  is  not  present, 
and  the  only  uncertainties  are  the  normal  and  malignant 
tissue  properties. 

Let  4  =  0,  from  Eq.  59  and  Eq.  61,  we  have 

<&  =  AeTAj(A1CA'[)-1A1Ae 
4  =  A6T(B1Ai)T(B1A1C(B1A1)T)~1B1A1Ae. 

(65) 

The  inequality  relationships  can  be  inherited  from  Eq.  64 
and  Eq.  62.  We  have 


dl>4n>  dl  (66) 

2.  Case  2:  Assume  signal  er  under  binary  hypotheses  are 
determined,  in  stead  of  random.  The  uncertainty  in  this 
problem  is  the  additive  measurement  noise  only. 

Let  C  be  Zero  matrix  and  let  A  replace  Ai,  from  Eq.  59 
and  Eq.  61,  we  have 

4n  =  AerAe/<T2 

4  =  AeTBf(B1Bf)-1B1Ae/4  (67) 

A©  =  A  AO. 

The  inequality  relationship  can  be  inherited  from  Eq.  64, 
and  we  have 

4n>d2r.  (68) 

The  forms  of  the  deflection  coefficients  and  their  inequality 
relationships  for  the  two  special  cases  can  also  be  derived 
directly  from  the  corresponding  special  binary  hypotheses. 

Discussion 

The  form  of  the  deflection  coefficients  for  special  case  1, 
and  their  inequality  relationships  demonstrate  that  it  is  the 
reduction  of  the  random  variable  space  during  the  forward 
and  reconstruction  procedures  that  causes  the  decrease  of 
the  detection  performance.  The  detection  performance  of 
the  optimal  detector  using  the  reconstruction  data  can  not 
be  better  than  that  using  the  measurement  data,  even  with¬ 
out  any  measurement  noise. 

The  form  of  the  deflection  coefficients  for  special  case  2 
demonstrates  that  dm  and  dr  are  inversely  proportional  to 
cr£,  which  reflects  the  direct  effect  of  the  measurement  noise 
on  the  detection  performance.  In  addition,  the  forms  and 
their  inequality  relationships  also  demonstrate  that  the  in¬ 
adequate  rank  of  the  reconstruction  matrix  B,  which  is  less 
than  the  length  the  measurement,  data,  will  reduce  the  ran¬ 
dom  variable  space  in  dealing  with  the  uncertainty  of  the 
measurement  noise,  and  cause  the  decrease  of  the  detec¬ 
tion  performance  at  the  reconstruction  domain.  The  effect 
of  the  measurement  error  on  the  detection  performance  has 
been  negatively  strengthened  by  the  reconstruction  proce¬ 
dure.  This  conclusion  agrees  well  with  the  simulation  result 
in  [16]. 

The  more  general  form  of  dm  and  dr  in  Eq.  59  and  Eq. 
61  clarify  the  effects  of  the  combination  of  the  projection 
and  reconstruction  procedures  as  well  as  the  measurement 
noise.  The  inequalities  in  Equations  64  and  62  illustrate 


^the  degradation  of  the  detection  performance  during  the 
processing  procedures.  This  general  conclusion  is  consis¬ 
tent  with  our  simulation  results  for  the  performance  of  the 
microwave  imaging  system. 

Clarifying  the  aspects  that  causes  the  degradation  of  the 
detection  performance  could  help  us  design  a  better  sys¬ 
tem.  It  suggests  that  good  sampling  of  the  measurement 
data,  which  maximizes  the  random  variable  space  in  deal¬ 
ing  with  the  uncertainties(i.e.  increase  the  rank  of  A),  will 
improve  the  detection  performance.  If  the  measurement  er¬ 
ror  cannot  be  ignored,  the  length  of  the  measurement  data 
should  be  limited  (i.e.  decrease  the  difference  in  rank  of  B 
and  Im)  to  reduce  the  degradation  of  post-processing  the 
reconstructed  image.  It  also  suggests  that  making  optimal 
decisions  before  reconstruction  will  improve  the  detection 
performance. 
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ADDENDUM  to  Final  Report  for  Award  Number  DAMD1 7-01 -1-05 13 
Principal  Investigator:  Loren  W.  Nolte,  PhD 

Contract  Organization:  Duke  University,  Durham,  North  Carolina  27708-0291 
Date:  December  17,  2003 

l.(a).  In  the  final  report  the  following  items  should  be  added  to  the  “Reportable 
Outcomes” 

Liewei  Sha,  Loren  W.  Nolte,  Zhong  Qing  Zhang,  and  Qing  Huo  Liu,  “Incorporating 
statistical  models  of  tissue  permittivity  into  an  optimal  signal  detection  theory  framework 
for  the  detection  of  breast  cancer  at  microwave  frequencies,”  submitted  to  the  IEEE 
Transactions  on  Medical  Imaging,  December,  2003.  (Copy  of  submitted  full  paper 
attached.  Note:  Funding  provided  by  the  USAMRMC  contract  has  been  acknowledged 
in  this  full  paper.) 

This  paper  should  also  be  added  to  the  “Bibliography  of  all  publications”  of  the  final 
report  as  item  5. 

(b) .Liewei  Sha,  a  PhD  student  supported  in  part  by  this  award,  expects  to  receive  her  PhD 
in  Spring,  2004. 

(c) .  Liewei  Sha  received  an  offer  of  employment  from  GE  Medical  Systems  in 
Milwaukee,  Wisconsin,  during  the  second  year  of  this  award,  based  on  experience 
supported  by  this  award. 
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2.  Research  highlights  completed  in  the  second  year. 

During  the  second  year  a  more  in  depth  treatment  of  the  signal  detection  theory  approach 
was  developed,  a  broader  class  of  possible  tumor  boundaries  were  considered,  lower 
order  approximations  were  developed  for  the  receiver  operating  characteristic  (ROC)  and 
the  probability  of  correct  localization  (PCL)  to  decrease  computational  complexity, 
convergence  issues  of  Monte  Carlo  approaches  were  quantified,  deflection  coefficient 
(delectability  index)  inequalities  were  developed  to  gain  insight  into  tumor  detection 
performance  tradeoffs,  and  the  impact  of  sensor  configuration  on  tumor  detection 
performance  was  determined  using  the  ROC. 

More  specific  results  regarding  these  topics  are  contained  in  the  recently  submitted  full 
paper  given  above,  referred  to  in  this  paragraph  as  Ref.  5,  i.e.  Liewei  Sha,  Loren  W. 
Nolte,  Zhong  Qing  Zhang,  and  Qing  Huo  Liu,  “Incorporating  statistical  models  of  tissue 
permittivity  into  an  optimal  signal  detection  theory  framework  for  the  detection  of  breast 
cancer  at  microwave  frequencies,”  submitted  to  the  IEEE  Transactions  on  Medical 
Imaging,  December,  2003.  A  more  in  depth  development  of  the  signal  detection  theory 
approach  and  how  the  measurements  and  models  of  breast  tissue  uncertainty  are 
incorporated  is  presented  in  Ref.  5,  section  II  A,  B  and  section  III.  Incorporating  tumors 
with  different  possible  boundaries  was  incorporated  into  the  algorithms.  In  particular,  the 
ROC  and  the  PCL  performance  measures  were  compared  for  tumors  with  sharp 
boundaries  versus  smooth  (Ref.  5,  section  Y).  Lower  order  approximations  were 
developed  for  the  Bayesian  processors  that  incorporate  the  uncertain  permittivity  for  both 
detection  and  localization.  Quantitative  tumor  detection  and  localization  performance 
measures,  the  ROC  and  the  PCL,  were  obtained  for  these  lower  order  approximations,  the 
EPB  UPI  and  LPB  SEUPI.  That  work  enables  one  to  greatly  reduce  the  computational 
time  for  the  ROC’s  and  PCL  performance  curves  for  both  tumor  detection  and 
localization  (Ref  5,  section  IV).  Since  Monte  Carlo  integration  is  necessary  in  both  the 
Bayesian  algorithms  as  well  as  in  the  lower  order  approximations,  the  convergence  of 
these  algorithms  was  investigated  in  order  that  no  more  iterations  than  are  necessary  are 
used  to  get  good  approximate  results  of  tumor  detection  performance  (Ref  5,  section  Y). 
Since  most  of  the  quantitative  detection  performance  results  are  obtained  using  Monte 
Carlo  methods,  deflection  coefficient  (delectability  index)  inequalities  were  developed  to 
gain  analytical  insight  into  tumor  detection  performance  tradeoffs  (Ref.  5,  Appendix). 
The  impact  of  the  sensor  array  configuration  on  tumor  detection  performance  as  a 
function  of  relative  location  of  a  tumor  to  the  array  sensors  was  quantified  (Ref  5, 
section  F). 
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3.  Please  add  the  following  references  to  the  final  report.  (A  more  comprehensive  list  of 
references  is  contained  in  Reference  5  of  the  modified  final  report  “Bibliography  of  all 
publications”) 
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Mines  and  Minelike  Targets,  Orlando,  FL,  April  1999. 

A.  M.  Richardson  and  L.  W.  Nolte,  “A  posteriori  probability  source  localization  in  an 
uncertain  sound  speed,  deep  ocean  environment”,  Journal  of  the  Acoustical  Society  of 
America,  Vol.  89,  No.  5,  pp.  2280-2284,  May,  1991. 
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4.  Note:  Although  the  first  annual  report  period  covered  the  period  August  1,  2001  - 
July  31,  2002,  the  following  item  presented  during  the  second  year,  August  1,  2002  - 
August  30,  2003  of  the  no-cost  project  extension  was  listed  under  the  reportable 
outcomes,  references  and  appendices  sections  in  the  first  annual  report. 

Liewei  Sha  and  Loren  W.  Nolte,  “Computer-aided  algorithms  for  breast  tumor  diagnosis 
using  microwave  diffraction  measurements”.  Era  of  Hope  2002  Department  of  Defense 
Breast  Cancer  Research  Meeting,  Orlando,  Florida,  September  25-28,  2002.” 

In  the  first  annual  report,  this  should  have  been  listed  as  “a  poster  to  be  presented”  and 
the  “Copy  of  the  paper  by  Liewei  Sha  and  Loren  Nolte”,  item  3  in  the  appendix  of  the 
first  year  report  should  have  read  “Copy  of  the  poster 
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